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1.0  INTRODUCTION 


The  Stability  and  Control  Program  (SACP)  is  a  computerized 
set  of  methods  for  predicting  stability  and  control 
characteristics  of  conceptual  level  flight  vehicles.  The  program 
was  developed  as  an  in-house  effort,  by  the  design  branch 
( AFWAL/FIAD ) ,  to  fill  a  deficiency  in  the  stability  and  control 
analysis  area. 

The  methodology  used  in  SACP  is  a  combination  of  existing 
methodology,  modified  existing  methodology,  and  newly  created 
methods.  The  program  is  divided  into  three  major  modules:  1: 

aerodynamics  and  longitudinal  coef f icients ,  2)  stability 
coefficients  and  lateral  coefficients,  and  3)  stability  analysis. 
The  program  is  set  up  to  run  on  the  design  branches  PRIME  750 
computer  and  uses  TEKTRONIX' s  PLOTIO  IGL  graphics  software  for 
graphical  output. 


Volume  I  of  this  report  documents  the  methods  used  by  the 
program.  Volume  II  is  a  user's  manual  and  computer  program 
description  for  SACP. 


2.0  AERODYNAMIC  METHODS 


The  methods  described  here  compute  the  critical  Mach  number, 
aerodynamic  center,  and  the  longitudinal  stability  derivatives  for 
the  specified  aircraft.  The  coefficients  computed  include  lift, 
drag,  pitching  moment,  the  variations  in  these  with  angle  of 
attack,  and  the  dynamic  pressure  ratios.  Methods  are  used  to 
compute  coefficients  for  untrimmed,  trimmed,  and  high  lift 
conditions  ensuring  a  good  coverage  of  aircraft  configurations. 
Each  method,  and  its  use,  is  discussed  below. 

2.1  CRITICAL  MACH  NUMBER  ( MCR ) 

The  drag-divergence  Mach  number  or  Mach  critical,  is  defined 
as  the  Mach  number  at  which  a  rapid  drag  rise  interrupts  the 
subsonic  trend  in  drag.  The  British  method  of  predicting  the 
critical  Mach  number  for  two-dimensional  airfoils  as  described  by 
the  Royal  Aerodynamic  Society  (3)  appears  to  be  the  most  accurate 
empirical  method  available.  This  method  uses  the  Sinnott  "crest 
criteria",  where  the  low  speed  pressure  at  the  airfoil  crest  is 
related  to  the  drag-divergence. 

The  SACP  program  uses  a  method,  taken  from  the  Large  Aircraft 
Program  (1),  (2),  that  is  analogous  to  the  British  two- 
dimensional  critical  Mach  number  prediction  procedure  in  order  to 
predict  the  critical  Mach  number  for  a  finite-aspect-ratio  swept 
wing.  The  critical  Mach  number  is  defined  as  the  value  of  free 
stream  Mach  number  which  produces  a  local  supersonic  flow  measured 
normal  to  the  sweep  of  the  isobar  at  the  crest.  The  local  Mach 
number  normal  to  the  crest  isobar  is  taken  to  be  1.02  for 
conventional  airfoils  and  1.05  for  supercritical  airfoils  in  order 
to  define  free  stream  critical  Mach  number.  The  sonic  condition 
at  the  crest  can  be  predicted  by  means  of  a  simple  equation  in 
which  the  incompressible  pressure  at  the  crest  of  the  airfoil  and 
compressibility  factors  are  used. 

The  value  of  1.02  local  Mach  number  for  the  weak  shock  at 
crest  condition  for  drag  rise  used  in  Sinnotts  transonic  airfoil 
theory  (4)  was  established  empirically.  The  Royal  Aerodynamic 
Society  (3)  shows  that  this  method  should  predict  critical  Mach  to 
within  +  0.015  for  the  majority  of  conventional  two-dimension 
airfoils.  However,  as  shown  in  reference  (3),  with  "peaky" 
airfoils  (as  in  the  supercritical  airfoil)  the  onset  of  rapid  drag 
rise  may  be  delayed  until  the  shock  is  substantially  downstream  of 
the  crest.  The  predicted  value  of  critical  Mach  based  on  a  local 
Mach  of  1.02  at  the  crest  may  thus  be  conservative  by  more  than 
0.02,  and  a  local  Mach  of  1.05  is  needed  to  achieve  good 
correlation . 

The  following  sections  discuss  the  methods  used  to  predict 
the  pressure  distribution  around  an  airfoil  and  to  determine 
critical  Mach  from  the  pressure  at  the  crest  along  with  the  method 
used  to  estimate  the  drag  rise  above  critical  Mach. 
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2.1.1  Pressure  Coefficient  Calculations 

The  incompressible,  inviscid  pressure  distribution  around  the 
airfoil  is  first  defined.  From  this,  the  pressure  at  the  crest 
can  be  determined.  The  method  of  Weber  (5)  was  used  for  the 
pressure  distribution  calculations.  This  method  requires  the 
airfoil  surface  coordinates  to  be  determined  at  the  chordwise 
locations  defined  by 


XU(IV)  =  .5  *  (1  +  C0S(IV  *  PI/N)); 


(2.1.1) 


Where  IV  takes  on  the  integer  values  between  zero  and  N,  and 
N  may  be  any  integer,  in  this  program  N  is  set  equal  to  32. 

The  Weber  formula  is  essentially  a  second-order  linear  theory 
whereby  the  pressures  are  determined  from  multiplication  of  the 
matrix  of  thickness  and  camber  ordinates  of  the  airfoil  by  a 
matrix  of  constants  given  by  the  Royal  Aerodynamic  Society  (5). 

The  formula  for  the  incompressible  pressure  distribution  on 
an  infinite  sheared  wing  was  obtained  from  the  incompressible  form 
of  equation  93  in  reference  (5)  resulting  in 

1  -  CP  =  (1/(1  +  ( ( SA ( IV  )  +  SE( IV ) ) /COS  L)2)) 

*  ((COS  A  *  COS  L  *  (1  +  SA(IV)  +  SD(IV)) 

+  SIN  A  *  COS  L  *  (1  +  SC(IV)/COS  L)  *  SQRT( ( 1  -  IV)/IV))2 
-  (COS  A  *  SIN  L  *  (SA(IV)  +  SD(IV)) 

2 

+  SIN  A  *  SIN  L  *  (1  +  SC(IV)/C0S  L)  *  SQRT(  (  1  -  IV ) /IV )  )  ) 


+  SIN  L  *  COS 


2  A  *  (1  -  (1/(1  +  ( ( SB ( I V  )  +  SE( IV ) ) /COS  L )  ‘ 


(2.1.2) 


For  CP  on  the  upper  surface  the  +  is  used,  and  for  the  lower 
surface  CP  the  -  is  used.  Also, 


SA(IV)  =  ,  S1(IU,IV)  *  ZT(IU) 

IU=  1 


•  >■ „■  ».  v  v  V  • 


-  '.**3 


SB( IV) 


S2 ( IU, IV )  *  ZT(IU) 


SC(IV)  =  S3 ( IV, IV )  *  ZT(IU)  +  S3 ( NX, IV )  *  SQRT(R0C/2) 


N 

iuii 

N 

IU=1 

N 

IU=1 

N 

IU=1 


SD(IV)  =  .  S4  ( IU ,  IV )  *  ZS(IU) 


SE(IV)  =  y ;  S5 ( IU, IV )  *  ZS(IU) 


Where  the  SI,  S2,  S3,  S4,  and  S5  matrices  are  constants  given  by 
the  Royal  Aerodynamic  Society  (5).  ZT(IU)  and  ZS(IU)  are  the 
thickness  and  camber  distribution  at  the  control  point  IU  given  by 


zt(iu)  =  . 5 ( yu  +  yl; 

and 
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ZS(IU)  =  . 5 ( YU  -  yl: 


Where  YU  and  YL  are  the  upper-  and  lower-surface  ordinates  at  the 
control  point  IU  given  by 

X(IU)  =  .5(1  +  COS(IU  *  PI/N)) 

2.1.2  Critical  Mach  Number  Calculation  from  Crestline  Pressure. 


The  procedure  followed  to  determine  critical  Mach  for  swept 
wings  is  similar  to  the  procedure  outline  in  reference  (3)  to 
predict  the  critical  Mach  for  two-dimensional  sections.  In  the 
procedure.  Equation  2.1.1  is  used  to  compute  the  pressure 
distribution  around  the  airfoil  at  a  sweep  angle  determined  from 


1 

the  actual  computations  for  ZT  and  ZS  are  dependent  on  the 
type  of  airfoil  represented 
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L  =  ARCOS  (COS(C2SWEEP)  )n;  WHEN  C2SWEEP  <_  40 

L  =  ARCOS  ( COS ( C2SWEEP ) )n  +  .76604n)/2;  WHEN  C2SWEEP  >  40 

(2.1.3) 

where  C2SWEEP  is  the  mid  chord  sweep  at  the  semi-span  of  the  wing 
and  the  factor  n  is  determined  from 

n  =  ARW/ (1.4+  ARW ) 

and  ARW  is  the  aspect  ratio  of  the  wing. 


The  angle  L  represents  an  effective  isobar  sweep  angle  at  the 
mid  span  region  of  the  wing.  The  procedure  used  to  determine 
critical  Mach  based  on  the  crest  pressures  is  as  follows: 

1.  Determine  a  chordwise  incompressible,  inviscid  pressure 
distribution  for  an  angle  of  incidence  (A).  Integrate  the 
pressure  distribution  to  obtain  the  lift  coefficient  (CLI) 

2.  Determine  the  chordwise  position  of  the  crest  for  each  A,  the 
crest  being  defined  as  the  point  at  which  the  airfoil  surface  is 
tangential  to  the  undisturbed  freestream  direction. 

3.  Determine  the  incompressible  pressure  coefficient  at  the 
crest  ( CPCRST  )  . 

4.  Use  CPCRST  to  determine  critical  Mach  from  the  relation 

(P/PT)  *  (1  +  .02  *  MCR2  »  COS2  L)3,5  -  1 

CPCRST  =  0.7  *  MCR2/SQRT(1  -  MCR2  *  COS2  L)  (2.1.4) 

where  (P/PT)  is  the  ratio  of  local  static  pressure  to 
freestream  stagnation  pressure  as  determined  by 

P/PT  =  ( 1  +  0.2  *  ML2  )3- 5 

where  ML  is  the  local  Mach  normal  to  the  isobar  sweep  L  at 
the  crest  of  the  airfoil.  ML  is  set  equal  to  1.02  for  conventional 
airfoils  and  1.05  for  "peaky"  or  supercritical  airfoils.  Equation 
2.1.4  uses  a  Prandt 1-Glauert  compressibility  factor  to  correct  the 
incompressible  pressure  coefficient  for  Mach  number  rather  than 
the  Karman-Tsien  factor  used  in  reference  3.  Miller  (6)  and  the 
Royal  Aerodynamic  Society  (7)  recommend  using  the  Prandtl-Glauert 
factor  instead  of  the  Karman-Tsien  factor  in  the  critical  Mach 
prediction  method  for  highly  cambered  airfoils  or  general  airfoils 
at  high-lift  coefficients.  The  relationship  determined  by 
Equation  2.1.4  is  plotted  in  Figure  2.1.1. 


5.  Use  the  Prandtl-Glauert  compressibility  factor  ( BMDD2  ) 
evaluated  at  critical  Mach  to  obtain  the  lift  coefficient 
CLMDD  from 

CLMDD  =  CLI/BMDD2 

6-  Repeat  Steps  1  through  5  for  a  set  of  incidence  in  order  to 
obtain  a  drag-rise  boundary  from  the  set  of  points  (CLMDD,  MCR )  . 


The  critical  Mach  number  predicted  by  the  above  six  steps  is 
prevented  from  exceeding  the  critical  Mach  number  of  the  fuselage 
alone  (shown  plotted  in  figure  2.1.2).  For  aircraft  that  are  not 
area-ruled,  where  the  isobars  are  allowed  to  upsweep  at  the  wing- 
fuselage  juncture,  the  method  would  tend  to  over  predict  MCR  when 
the  value  approaches  the  fuselage  MCR.  The  prediction  -  versus  - 
test  MCR  correlation  shown  in  figure  2.1.3  is  thus  applied  for 
conventional-wing  predictions. 
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2.2  AERODYNAMIC  CENTER 

The  methods  for  calculation  of  aerodynamic  center  locations 
are  those  of  the  DATCOM  (8).  They  have  largely  been  taken  from 
the  Large  Aircraft  Program  ( 1 )  with  modifications  made  to 
incorporate  more  complex  airframes. 

2.2.1  Aerodynamic  Center  of  Forebody 

The  subsonic  location  of  the  aerodynamic  center  for 
forebodies  with  ogive  nosecones  is  approximated  in  the  DATCOM  (8) 
as 


XACNO 


-0.54  *  (BLN  +  1.6 


( XCREW  -  BLN) ) /CREW 


(2.2.1  ) 


where 

CREW  is  the  wing  exposed  root  chord  length 
BLN  is  the  body  nose  length 

XCREW  is  the  x  location  of  wing  exposed  root  chord  leading  edge 

Figure  2.2.1  defines  the  geometric  parameters  BLN,  XCREW,  and 
CREW.  The  supersonic  forebody  aerodynamic  center  is  obtained  from 

XACNO  =  XCREW/ CREW  *  ( XCPOL  -  1)  (2.2.2) 

where  the  term  XCPOL  is  obtained  from  figure  2.2.2 

2.2.2  Aerodynamic  Center  of  Surface  (Single  Panel). 

The  aerodynamic  center  of  exposed  surfaces  is  determined  from 
the  DATCOM  charts  presented  in  figures  2.2.3a  through  2.2.3f. 

These  charts  are  valid  for  subsonic  Mach  numbers  less  than  Mach 
critical  and  supersonic  Mach  numbers  greater  than  1.2.  For 
transonic  conditions,  the  data  presented  in  the  DATCOM  in  terms  of 
transonic  similarity  parameters  (Figures  2.2.4a  through  2 . 2 . 4d  ) 
are  used  to  determine  the  aerodynamic  center  position. 

The  procedure  for  obtaining  aerodynamic  center  can  be 
summarized  as: 


For  Mach  <_  MCR  or  Mach  _>  1.2 
XACR  =  XAC 1 


SPAN/2 


I*  Jl*, 


.6 


‘A  -< 


XCPOL 


.5 


fN  H 


.4 


(fraction  of  j. - - 

body  length) 

2 - 

.1 - 

o4— L-4 

0  .2 


Figure  2.2.2  Super 
witrh 


XAC1 

/x.  _\  -8* 


ARcanA^E  ; 

-  I  ■  -i  ■  — i- 

nrrn 


-6  *  ~  3  ~ 

=*7± 


(c)  X  =0.25  i 


III' 


i  1  :  I  L 

~  UNSWEPT  T.E 


ARtan.V 


a 


!  _L  i 

I  I  i  ii _ j_ 

i  !  I  SUBSONIC 


0  TAN  Ale  1 


_ 3_ 

TAN  A, 


!  I  I  I 

—  SUPERSONIC 


1  TAN  A, 


TAN  A, 


r 

(d)‘  X 

=0.33 

.  1  ! 

ARtanAi 

j 

LE  ! 

'  6 

Is! 

_ 

"-Li  !  - 

■ — 

_A  '  '  - - 

3  i  '  i 

j 

L 

i  i  ; 

i  i  i  i 

! 

1  • 

i 

MM 

! 

TAN  A, 


i  i  !  | 

UNSWEPT  T.E. 


SUBSONIC 

_ 0_ 

TAN  A. 


1  jSm 

IIms 


X  ARtanA 

j  j  !  _ 

1  ;  a 

nr 


i 


1  i  i 
1  !  - 

-  SUPERSONIC 

■  .  t 

I 

3 

TAN  A.  c 


TAN  A 


Figure  2.2.3  (Cont  *d) 


For  MCR  +  .05  >  Mach  >  MCR 


XACR  =  XAC 1  +  ( XAC2  -  XAC1 )  *  (Mach  -  MCR)/. 05 
For  FM2  >  Mach  >  MCR  +  .05 
XACR  =  XAC2 


For  1.2  >  Mach  _>  FM2 

XACR  =  XAC 1  +  ( XAC2  -  XAC1  )  *  (Mach  -  FM2)/(1.2  -  FM2  )  (2.2.3) 

where 

FM2  SQRT  (1  +  TOCS’6667) 

MCR  is  the  critical  Mach 

TOCS  is  the  thickness  to  chord  ratio  of  the  surface 

XAC 1  is  given  by  figure  2.2.3 

XAC2  is  given  by  figure  2.2.4 

2.2.3  Aerodynamic  Center  of  Multi  Panel  Surfaces. 

The  prediction  of  the  aerodynamic  center  location  of  multi 
panel  surfaces  is  taken  from  the  method  developed  at  Convair 
Aerospace  as  reported  by  Henderson  (9).  The  non-straight-  tapered 
wing  is  divided  into  panels,  with  each  panel  having  conventional, 
straight-tapered  geometry.  The  individual  lift-  curve  slope  and 
aerodynamic  center  are  estimated  for  each  panel,  using  the 
technique  described  above  for  the  single  panel  wing  and  treating 
each  constructed  panel  as  an  isolated  wing.  The  individual  lift 
and  aerodynamic  center  locations  for  each  constructed  panel  are 
then  combined,  using  an  inboard-outboard  weighted-area 
relationship 

XASCR  = 

XACR( 1  )*CLALPSO( 1  )*SXX0( 1  )  +  XACR ( 2 ) *CLALPS0 ( 2 ) »SXX0( 2 )  +... 
CLALPS0(  1  )  *  SXXO(I)  +  CLAI.PSO  ( 2  )  *  SXX0(2)  +  ... 

(2.2. 

where  the  outboard  wing  aerodynamic  center  is  referenced  to  the 
total  surface  root  chord  length  given  by 


XACR  =  ((XACR  *  CRO)  +  ( XCRO  -  XCRES))/CRES 
and 


CLALPSO ( i  ) 
SXX0( i ) 

CRO 

XCRO 


is  the  zero  AoA  lift  curve  slope  of  the  panel 

is  the  area  of  the  panel 

is  the  panel  exposed  root  chord 

is  the  X  location  of  the  panel  root  chord 

leading  edge 


XCRES 

CRES 


is  the  X  location  of  the  surface  root  chord 
leading  edge 

is  the  surface  exposed  root  chord 


The  geometry  for  the  multi  panel  arrangement  is  illustrated  in 
figure  2.2.5. 


2.2.4  Aerodynamic  Center  of  Wing-Lift  Carryover  on  the  Bod\ 


The  location  of  the  aerodynamic  center  due  to  the  wing-lift 


1 


carryover  on  the  body  is  determined  by  use  of  the  DATCOM  method. 
For  BARE  _>  4  the  subsonic  aerodynamic  center  location  is  obtained 
from 


XACBW  =  .25  +  DXQC  *  FDOB  *  TAN(C4SWEEP) 


DXQC  =  ( SPANW  -  BDW ) / ( 2  *  CREW )  (2.2.5) 

where 

SPANW  is  the  span  of  the  wing 

BDW  is  the  body  diameter  along  the  wing  root  chord 

CREW  is  the  wing  exposed  root  chord  length 

C4SWEEP  is  the  wing  quarter  chord  sweep 

BARE  is  the  Prandt 1-Glauert  compressibility  factor 

times  the  effective  aspect  ratio 

and  FDOB  is  determined  from  the  plot  in  figure  2.2.6. 

For  BARE  <  4  the  aerodynamic  center  location  is  determined  from 

XACBW  =  ( XACBWO  -  XACBW)  *  (BARE  -  4  ) 2/ 1 6  +  XACBW  (2.2.6) 

where  XACBW  is  from  the  equation  (2.2.5)  and  XACBWO  is  determined 

t>y 


XACBWO  =  .5  *  .25  *  ARWE  *  TAN ( LESWEEP )  *  (l  +  TRWE )  (2.2.7) 

and 

ARWE  is  the  exposed  aspect  ratio  of  the  wing 

TRWE  is  the  exposed  taper  ratio  of  the  wing 

LESWEEP  is  the  leading  edge  sweep  of  the  wing 
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Figure  2.2.5  Geometry  for  the  Multi-Panel  Surfaces 
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Equation  2.2.7  is  limited  to  values  less  than  or  equal  to 
0.5.  For  supersonic  conditions,  the  aerodynamic  center  location 
is  estimated  from  figure  2.2.7.  For  transonic  conditions,  the 
aerodynamic  center  location  is  determined  by  linear  interpolation 
of  the  aerodynamic  center  values  determined  at  the  critical  Mach 
number  and  Mach  1.1. 


2.2.5  Total  Aerodynamic  Center. 


With  the  above  computed  values,  the  location  of  aerodynamic 
center  for  a  wing  body  configuration  is  given  by  the  DATCOM  as 


XACWB  +  KS  *  XACS(i) 


XAC  = 


1 


KS 


(2.2.8 


where  the  XACS(i)  are  the  aerodynamic  centers  of  the  aircraft 
components 


CLAS 


KS 


CLAWB 


NS 

NW 


AS 

*  AW  * 


COS(DIHS) 


( 1  -  DEWS 


(2.2.9) 


and 


CLAS 

CLAWB 

NS 


NW 


AS 

AW 

DIHS 


is  the 
is  the 
is  the 
section 
is  the 
section 
is  the 
is  the 
is  the 
is  the 
wing  ( f 


zero  alpha  lift  curve  slope  for  the  surface 
zero  alpha  lift  curve  slope  for  the  wing-body 
dynamic  pressure  ratio  at  the  surface  (from 
2.8  ) 

dynamic  pressure  ratio  at  the  wing  ( from 

2.8  ) 

exposed  area  of  the  surface 

exposed  area  of  the  wing 

dihedral  angle  of  the  surface 

downwash  factor  for  the  surface  from  the 

rom  section  2.6) 


Equation  2.2.9  is  summed  over  all  of  the  surfaces  on  the  airframe. 
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Ref.  8  ,  FIGURE  4.3.2.1-37 

(»)  WITH  AFTERBODY 


EFFECTIVE  RATIO  OF  BODY  DIAMETER  TO  ROOT  CHORD,  — 


Figure  2.2.7  Aerodynamic-Center  Locations  for  Lift  Carry¬ 
over  onto  Body  at  Supersonic  Speeds 
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2.3  ZERO-LIFT  PITCHING  MOMENT  ( CMO ) 


The  method  for  predicting  the  zero-lift  pitching  moment  for 
an  aircraft  configuration  considers  only  the  effects  of  surface 
components  on  CMO.  It  does  not  include  the  effect  of  an 
asymmetrical  fuselage  or  the  effect  of  stores  and  nacelles  located 
near  the  surfaces. 

The  subsonic  zero-lift  pitching  moment  for  lifting  surfaces 
with  linear  twist,  up  to  the  critical  Mach  number,  is  given  in  the 
DATCOM  as 


CMOS  =  ( CMON  +  CMOOT  *  TWIST)  *  CMACH 


where  CMON  is  the  CMO  of  an  untwisted  surface  given  by 


(2.3.1) 


CMON  =  (ARS  *  COS2 ( C4 SWEEP ) ) / ( ARS  +  2  *  COS ( C4SWEEP ) )  *  CMOSECT 

(2.3.2) 

CMOSECT  is  the  average  section  pitching  moment  coefficient 

determined  by  averaging  the  section  CMO  for  each 
surface  panel  using 


CMOSECT 


CMOOT 


CMACH 


TWIST 

ARS 

C4SWEEP 

SCAMP 

SXXS 

CMC4 

ASEXS 

TOCS 

MACH 


( SCAMP ( i )  *  SXXS(i)  *  CMC4 ( i  )  ) / ASEXS 


(2.3.3) 


is  determined  by  a  table  lookup  from  lifting  line 
theory  graphs  as  shown  in  figure  2.3.1 

(1  +  5.9  *  TOCS  *  MACH5 )/(SQRT( 1  -  MACH2  COS2 ( C4SWEEP ) ) 

(2.3.4  ) 


is  the  twist  on  the  surface 

is  the  aspect  ratio  of  the  surface 

is  the  quarter  chord  sweep  of  the  surface 

is  the  camber  on  the  panel 

is  the  exposed  area  of  the  panel 

is  a  constant  based  on  the  airfoil  type 

is  the  exposed  area  of  the  surface 

is  the  thickness  to  chord  ratio 

is  the  current  Mach  number  (not  greater  than  Mach 
critical ) 


The  CMO  for  the  total  configuration  can  then  be  computed  by 
summing  the  individual  CMOs  for  surfaces  after  they  have  been 
converted  to  the  reference  area.  The  equation  is 


CMO  =  '  .CMOS ( i  )  *  S ( i  )  / SERF 


(2.3.5) 


where 

S(i) 

SREF 


is  the  area  of  the  surface 
is  the  reference  area 
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2.4  ZERO  ALPHA  LIFT  COEFFICIENT  ( CLO ) 


The  method  for  predicting  the  zero  alpha  lift  coefficient  has 
been  taken  from  the  Large  Aircraft  Program  (LACP)  (1).  the  method 
first  computes  the  lift  curve  slope  and  the  ALPHA  where  zero  lift 
occurs.  CLO  is  then  computed  as  the  product  of  these  two  values; 

CLO  =  CLA  *  -ALO 

This  equation  holds  for  all  surfaces,  bodies,  and  nacelles. 

2.4.1  Lifting  Surface  CLO 

For  each  surface  of  interest  the  lift  curve  slope,  CLA,  is 
defined  by 

CLA  =  CLAB  *  XT  (2.4.1) 

where 

CLAB  is  the  surface  alone  CLA  with  no  thickness  effects.  The 

factor  XT  accounts  for  the  effect  of  airfoil  thickness  plus 
camber.  The  equation  for  CLAB  was  originally  taken  from  Polhamus 
(10)  and  adopted  for  the  LACP.  The  final  form  is; 

CLAB  =  1.0/ (57. 3  *  ( ( M*/MACH ) Z/CLAO  +  BETAP/4.0));  WHEN  MACH  >  M* 

(2.4.2) 


and 


CLAB 


(.0548311  *  ARS)/[1.0  +  SQRTd.O  +  (1.0  -  COS( C2SWEEP ) 1 • 334 
*  (MACH/M*)2'667]  *  ( ARS/ (2.0  *  COS ( C2SWEEP ) ) ) 2 ) ) ; 

WHEN  MACH  <  M* 


(2.4.3) 


Where  M*,  a  function  of  ARS  and  C2SWEEP,  is  the  limiting  Mach, and 
is  defined  by; 

M*  =  CSUBO  +  (1.0  -  CSUBO )  *  (1.0  -  COS ( C2SWEEP ) ) 2 


and 


CSUBO 
BETAP 
Y  =  (  1 

Z  =  M* 


=  (10.0  +  0.91  *  ARS3  )/( 10.0  +  ARS3) 

=  (MACH  -  M* )  *  (1.0  +  (M*/MACH)Y)2 

+  PI  *  ARS )  / ( 3  +  PI  *  ARS)  *  (2  +  .667 

*  CLAO  +  ARS2/ (3  *  PI  *  ARS/CLAO  *  (PI 
COS ( C2SWEEP  )  • 667  ) 


*  SQRT ( TRS 

*  ARS/CLAO 


TRS2 
.  0  ) 


where 


CLAO 

C2SWEEP 

ARS 

TRS 


is  defined  as  CLAB  at  MACH  =  M* 

is  the  half  chord  sweep  of  the  theoretical  surface 
is  the  aspect  ratio  of  the  theoretical  surface 
is  the  taper  ratio  of  the  theoretical  surface 


‘  Wings  having  thijk  airfoils  undergo  a  degradation  in  CLA 

b  beginning  at  MACH  >  M  .  The  level  of  CLA  versus  Mach  dips, 

usually  reaching  a  minimum  at  MACH  <  1.0,  and  then  recovers  to  a 
second  peak  at  M  >  1.0.  To  account  for  this  phenomenon,  the  basic 
CLA  equations  have  been  modified  by  a  factor  KT,  as  defined  by; 


KT  = 

1.0- 

(4.0 

*  SIG1  * 

(1.0 

-  SIG1  )  )3 

*  GAMMA; 

Ml 

£ 

MACH 

< 

M2 

KT  = 

1.0- 

(4.0 

*  SIG2  * 

(1.0 

-  SIG2  )  )3 

*  GAMMA; 

M2 

MACH 

<_ 

M3 

KT  = 

1.0; 

Ml 

> 

MACH 

> 

M3 

(2.4.4) 

where  the  factors  GAMMA,  SIG1 ,  SIG2,  Ml,  M2  and  M3  are  defined  by 
the  equations  that  follow: 


GAMMA  =  9.0  *  DTOCL/ (1.0  +  .5  *  ARDT ) 

DTOCL  =  (TOCS  -  TOCL ) /COS ( C2SWEEP ) 

TOCL  =  1.0/ (4. 4  *  ARS  *  COS ( C2SWEEP ) 1  *  5  ) 

ARDT  =  ARS  *  DTOCL  (2.4.5) 

SIG1  =  .5  *  (MACH  -  M 1 ) / ( M2  -  Ml) 

SIG2  =  .5  *  (1.0  +  (MACH  -  M2)/(M3  -  M2) 

Ml  =  1.0  -  2.0  *  TOCS  *  ( ARS3 / (4.0  +  ARS3)) 

*  COS ( C 2 SWEEP )  ’ 5  *  (1.0  +  1.5  *  CLD1'5) 

M2  =  Ml  +  TOCS 

M3  =  1.0  +  TOCS  (2.4.6) 

and 

TOCS  is  the  thickness  to  chord  ratio  of  the  surface 

CLD  is  the  camber  of  the  surface 


v  jN 


NOTE : 


If  DTOCL  >  .07  then  DTOCL  =  .07 
If  ARDT  >  .1  then  ARDT  =  .1 


0  <_  Ml  <_  M 

If  Ml  >  M*  then  Ml  =  M* 


The  derivation  of  equation  2.4.4  is  based  on  the  data  trends 
and  analysis  of  Polhumas  (11)  and  on  other  limited  data  (e.g.. 
Nelson  and  McDevitt  (12),  and  Tinling  and  Kolk  (13)).  It  should 
be  noted  that  a  CLA  'bucket"  is  predicted  only  if  the  wing 
streamwise  airfoil  TOCS  exceeds  the  limit  thickness  defined  by 
equation  2.4.5.  The  limit  thickness  boundary  was  established  from 
the  statistical  boundaries  presented  by  Donlan  and  Weil  (14). 


The  angle  of  attack  at  zero  lift,  ALO  is  determined  by 


ALO  =  AALOT  *  TWIST 


(2.4.7  ) 


where 

TWIST 

AALOT 


is  the  twist  on  the  surface  (-  for  washout) 
is  obtained  from  a  curve  fit  of  parametric  data 
reported  by  Gilman  and  Burdges  (15)  for  wings  with 
linear  element  twist.  The  derived  equation  for 
AALOT  is; 


AALOT  =  .093  -  (.000571  *  SPBETA )  +  (.5761  *  TRW) 
-  (  .2645  *  TRS2 ) 


(2.4.8) 


where 

SPBETA 

BETA 


C4SWEEP 


ATAN( TAN( ABS ( C4 SWEEP /BETA ) ) ) 
SQRT( ABS ( 1  -  MACH2 ) ) 


is  the  quarter  chord  sweep  of  the  surface 


2.4.2  Body  and  Nacelle  CLO 


The  equations  for  calculation  of  CLA  for  aircraft  bodies  and 
nacelles  in  LACP  have  been  simplified  for  SACP.  CLA  is  now 
defined  as; 


CLA  =  2.0  *  CSA/SREF 


(2.4.9) 


where 

CSA 

SREF 


is  the  maximum  cross  sectional  area  of  the  component 
is  the  aircraft  reference  area. 
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ALO  in  this  case  is  also  simplified.  For  nacelle  components  ALO 
is  just  the  incidence  of  the  nacelle.  For  bodies  ALO  is  defined 
as 

ALO  =  AT AN ( ( BTL  *  TAN( BTUSA ))/( BL/2 . 0 ) )  (2.4.10) 

where 

BTL  is  the  boattail  lepgtfc 

BTUSA  is  the  boattail  upsweep  angle 

BL  is  the  length  of  the  body 

2.4.3  Total  Aircraft  CLP 

Once  all  of  the  component  CLO  terms  have  been  computed  and 
referenced  to  the  aircraft  reference  area,  then  the  total  CLO  is 
expressed  as  the  sum  of  all  of  the  components  as; 


CLOT  =  CLOB  +  CLON  +  CLOW  * 
+  CLOH  *  (SH/SREF)  + 


CLOB 

CLON 

CLOW 

CLOC 

CLOH 

CLOV 

CLOF 


is 

is 

is 

is 

is 

is 

is 


the 

the 

the 

the 

the 

the 

the 


CLO 

CLO 

CLO 

CLO 

CLO 

CLO 

CLO 


for 

for 

for 

for 

for 

for 

for 


the 

the 

the 

the 

the 

the 

the 


(SW/SREF)  +  CLOC 
CLOV  *  (SV/SREF) 

body 

nacelles 

wing 

canard 

horizontal  tail 
vertical  tail 
fins 


(SC/SREF) 

CLOF  *  (SF/SREF) 
(2.4.11 


The  effect  of  wing  incidence  is  not  included  in  the  CLO 
calculation.  The  wing  incidence  is  added  into  the  total  effective 
angle  of  attack  (See  section  2.6.1) 
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2.5  ZERO  LIFT  DRAG  ( CDO ) 


The  drag  of  an  aircraft  can  be  represented  as  the  sum  of 
minimum  drag,  plus  drag  due  to  lift,  plus  drag  due  to  trim.  The 
SACP  program  maintains  separately,  the  minimum  drag  (or  zero  lift 
drag)  and  the  trim  drag.  CDO  is  comprised  of  the  drag  items  that 
are  assumed  to  be  independent  of  lift,  such  as  friction,  form, 
interference,  wave,  base,  and  camber.  These  components  and  drag 
rise  are  discussed  below. 

A  large  part  of  the  subsonic  CDO  is  comprised  of  the  sum  of 
the  above  mentioned  items.  The  drag  of  each  component,  then,  is 
computed  as 

CD  =  CF  *  (AWET/SREF)  *  FF  *  IF  *  CDW  *  CDC 
for  surfaces  and 

CD  =  CF  *  (AWET/SREF)  *  FF  *  IF  *  CDW  *  CDS  (2.5.1) 

for  bodies  and  nacelles,  where 

CF  is  the  compressible  flat  plate  skin-friction 

coefficient 

AWET  is  the  component  wetted  area 

SREF  is  the  reference  area 

FF  is  the  form  factor 

IF  is  the  interference  factor 

CDW  is  the  wave  drag 

CDC  is  the  drag  due  to  camber  and  twist  of  surfaces 

CDB  is  the  base  drag  for  bodies  and  nacelles 

2.5.1  Friction  Drag  (CF) 

The  flat-plate,  compressible,  turbulent,  skin-friction 
coefficient  is  determined  from  methods  given  in  the  Large  Aircraft 
Program  ( 1 ) .  These  methods  are  based  on  work  by  White  and 
Christoph  (16)  and  give  the  following  expression  for  CF: 

CF  =  T  *  F2  *  . 430/ ( LOG , Q ( RNL  *  L  *  T1’67  *  F))2'56  (2.5.2) 

where 

T  =  1  . 0/ (  1  . 0  +  0.178  *  MACH2  ) 

F  =  1.0  +  0.03916  *  MACH2  *  T 

L  is  the  characteristic  length  of  the  component  (mean 

aerodynamic  chord  for  surfaces  or  length  for  bodies  and 
nacelles ) . 
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is  based  on  either  component  length  or  an  admissible 
surface  roughness,  whichever  produces  a  smaller 
value  of  Reynolds  number,  as  follows: 


RNL  =  MIN( ( RNF  *  L),  (K1  *  ( L/K ) 1  *  0489  )  ) 


(2.5.3) 


where 


is  the  Reynolds  number  per  foot  determined  from 
standard  atmospheric  tables 

is  admissible  surface  roughness  and  is  an  input 
quantity 


K 1  =  37.587  +  4.617  *  MACH  +  2.949  *  MACH2  +  4.132  *  MACH3 


For  mixed  laminar-turbulent  flow,  transition  location  is 
specified  for  the  upper  and  lower  surfaces  of  the  wing.  For  the 
laminar  portion  of  the  flow,  the  Blasius  skin-friction  relation  is 


CF(XR)  =  CFCFIL  *  1 . 328/SQRT( RNL  *  XR ) 


(2.5.4 


where 


CFCFIL 


=  (1.0  +  0.1256  *  MACH2)-*12,  is  used  up  to  the 
transition  point,  XR,  the  laminar  momentum  thickness, 
which  begins  some  fictitious  distance,  DX,  ahead  of 
transition.  The  skin-friction  coefficient  for  the 
turbulent  part  of  the  flow  is  calculated  from  equation 
(2.5.2),  where  the  Reynolds  number  is  calculated  from 


RNL  =  (DX  +-  L  -  XR)  *  RNF 

The  value  of  CF  with  transition  is  finally  given  by 
CF  =  (DX  +  L  -  XR)/L  *  CFTURB 


(2.5.5) 


(2.5.6) 


Calculated  values  of  CF  versus  RNL  are  presented  in  figures  2.5.1 
through  2.5.6  for  mixed  laminar-turbulent  flow. 

2.5.2  Form  Factors 

The  component  form  factors,  FF,  account  for  the  increased 
skin-friction  caused  by  the  super  velocities  of  the  flow  over  the 
body  or  surface  and  the  boundary-layer  separation  at  the  trailing 
edge.  The  form  factor  for  the  "body"  component  is  computed  as 


FF  =  1.0  +  60.0/FR3  +  0.0025  *  FR 


(2.5.7) 


where 
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FR  is  the  fineness  ratio  of  the  component. 

For  "nacelle"  components,  the  form  factor  is  given  by 

FF  =  1 . 0  +  0 . 35/FR  (2.5.8) 

Equations  2.5.7  and  2.5.8  were  obtained  from  the  Convair  Aerospace 
Handbook  (17)  and  also  appears  in  the  DATCOM  (8). 

The  airfoil  form  factors  depend  upon  airfoil  type  and 
streamwise  thickness  ratio.  For  6-series  airfoils,  the  form 
factor  is  given  by 

FF  =  1.0  +  1.44  *  TOCS  +  2.0  *  TOCS2  (2.5.9) 

where 

TOCS  is  the  surface  thickness  to  chord  ratio. 

For  4-digit  airfoils,  the  form  factor  is  given  by 

FF  =  1.0  +  1.68  *  TOCS  +  3.0  *  TOCS2  (2.5.10) 

For  biconvex  airfoils,  the  form  factor  is  given  by 

FF  =  1.0  +  1.2  *  TOCS  +  100.0  *  TOCS4  (2.5.11) 

And  for  supercritical  airfoils,  the  form  factor  is  given  by 

FF  =  1.0  +  ZK  *  (CLD/.4)  +  1.44  *  TOCS  +  2.0  *  TOCS2  (2.5.12) 

The  factor  ZK  *  CLD/ . 4  in  equation  2.5.12  is  an  empirical 
relationship  which  shifts  the  6-series  form  factor  equation  to 
account  for  the  increased  super  velocities  caused  by  the 
supercritical  section  design  camber  CLD.  The  factor  ZK/.4 
(derived  from  experimental  data)  is  shown  plotted  in  figure  2.5.7, 
as  a  function  of  the  Mach  number  relative  to  the  wing  Mach 
critical.  The  equations  above  were  extracted  from  the  Large 
Aircraft  Program  (l). 

2.5.3  Interference  Factors 


The  component  interference  factors,  IF,  account  for  the 
mutual  interference  between  components.  For  the  fuselage,  the 
interference  factor  is  given  by 

IF  =  RWB  ;  WHEN  MACH  <  1.0 

IF  =  RLS  *  RWB  ;  WHEN  MACH  _>  1.0  (2.5.13) 

where 
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RWB  is  shown  plotted  in  figure  2.5.8  as  a  function  of 

fuselage  Reynolds  number  and  Mach. 

For  other  bodies,  such  as  engine  nacelles,  the  interference 
factor  is  given  based  on  experimental  experience  with  similar 
configurations.  The  Convair  Aerospace  Handbook  (17)  recommends 
using 

IF  =  1.0  for  nacelles  mounted  out  of  the  local  velocity  field  of 
the  wing 

IF  =  1.3  for  nacelles  mounted  in  moderate  proximity  to  the  wing 

IF  =  1.5  for  nacelles  and  stores  mounted  flush  to  the  wing  or 
fuselage 

SACP  currently  requires  the  user  to  input  the  interference  factor 
for  nacelles. 

The  interference  factor  for  the  surfaces  is  computed  as 
IF  =  RLS 

where  RLS  is  the  lifting  surface  interference  factor  presented  in 
figure  2.5.9. 

If  the  Mach  number  is  greater  than  1.0  then  the  interference 
factors  for  both  bodies,  and  surfaces  are  set  to  1.0. 

2.5.4  Wave  Drag 

Supersonic  wave  drag  is  determined  on  the  basis  of  a 
component  buildup  for  which  simplified  shapes  are  assumed.  Three 
basic  simplified  shapes  are  used  to  represent  the  airplane: 
bodies,  nacelles,  and  surfaces.  The  component  buildup  assumes 
that  the  total  drag  is  the  sum  of  the  isolated  wave  drag  of  each 
component  and  does  not  allow  for  the  mutual  interference  between 
components.  However,  the  component  buildup  method  does  give  wave- 
drag  results  comparable  to  favorable  and  unfavorable  interference. 

2.5.4. 1  Surface  Wave  Drag 

The  technique  used  to  estimate  surface  wave  drag  evolved  from 
a  method  that  applies  transonic  similarity  theory  to  straight 
surfaces.  Data  correlations  at  Mach  1  were  performed  on  a  large 
number  of  unswept  surface  configurations  with  blunt  and  sharp 
leading-edge  airfoils.  For  SACP,  these  results  were  represented 
by  an  analytical  function  common  to  both  types.  The  equations, 
which  are  taken  from  the  Large  Aircraft  Program  (1 ),  were  modified 
for  M  >  1  to  produce  a  peak  value  at  low 


COS  Max  Thickness  :weep 


lifting  Surrace  Correlation  Fact:r  For 
Subsonic  Minimum  Drag 


supersonic  speeds  and  then  to  decrease  at  high  Mach  numbers  to 
values  predicted  by  straight-wing  linear  theory  for  equivalent 
two-dimensional  configurations.  Finally,  sweep  effects  were 
included.  The  resulting  semi-empirical  equations  are  presented 
below 


2  *  Kt  *  Kw  *  Kc  *  Kb  *  TOC 1,667 
BETAW  *  Kb  *  Kw  *  FBXM  +  T1  +  T2 


.1  .  667 


+  3 . 33  *  Kt  *  Kw  *  Kc  *  Kb  *  TOC'-oa' 
BETAW  *  Kb  *  KW1  *  FBXM  +  T3  +  T4 


2.5.14 


where 


is  airfoil  thickness  distribution  factor 


Kt  =  1.0  +  4.0  *(.5  -  XTOC  *(1.0  +  .5  *  SQRT(ROT)))2 
-  .25  *  SQRT(ROT)  *  (1.0  -  XTOC)2 


(2.5.15) 


Kb  and  Kw  are  airfoil  factors 


Kb  = 

Kb  = 
XTOC  - 
ROT  = 
Kc  = 


1.0,  Kw  =  1.2  for  double-wedge  sections 

1.069,  Kw  =  1.0  for  curved-type  sections 

is  the  location  of  airfoil  section  mean  Y  ordinate 

is  the  leading  edge  radius  to  thickness  ratio 

is  the  airfoil  camber  factor  given  by 


Kc  = 


1.0  +  2.5  *  HOT' 


CAMBER  *  . 055/TOC 


(2.5.16) 


Kp  =  ( COS ( LESWEEP  )  +  TR 1  * ( TAN ( TESWEEP  )  2  -  TAN ( LESWEEP 
(1  +  TR2  *  (TAN (LESWEEP)  +  TAN ( TESWEEP )) 2 ) 


)2  )  )/ 


(2.5.17 


TR  1  =  0 . 5 / (  1  . 0  +  ZLAM)2 

TR2  =  1 . 0/ ( 1 .0  +  ZLAM)2 
LESWEEP  is  the  leading  edge  sweep 
TESWEEP  is  the  trailing  edge  sweep 

ZLAM  is  the  exposed  taper  ratio 


T 1  =  (  1  . 
T2  =  ARE 


0/ARE)/(l.0  +(1.0  +  ZLAM)  *  FB  * 
3 / (  1.0  +  .33  *  ARE3  *  BETAW4) 


BETAW1 +Kw ] 


T3  =  ( 2/ARE3  )/(  1  . 0  +  (.  667  +  ZLAM)  *  FB  *  BETAW1+Kwl) 
T4  =  1.0/ (1.0  +  3.0  *  ARE  *  BETAW4) 


=  Kw- 


BETAW  =  SQRT(M^  -  l)/TOC-333 

FBXM  =  FBXM 

FB  =  .3  +  .7  *  Kp 

FEX  =  1.0  +  2.0  *  ZLAM’ 33 J 

XM  =  .5  *(1.0  +  ZLAM2  *(2.0  -  ZLAM)3) 


If  Beta  ( SQRT ( 


1  .  0  )  )  >  TAN( LESWEEP ) 


XM  =  (TAN( LESWEEP) /BETA)2 
Z  =  COS (LESWEEP)  +  COS(TESWEEP) 


where 

ARE  is  the  straight -surf ace  ARW  having  the  same  value  of 

CDW  at  M  =  1.0  as  ARW,  where; 

ARW  =  AR  * ( TOCS  )  *  333 

and 


AR 

M 

TOCS 


is  the  aspect  area  of  the  exposed  surface 
is  the  Mach  number 

is  the  thickness  ratio  of  the  surface 


The 

equation 


value  of  ARE  is  determined  by  solving  the  following 
using  an  iterative  method 


2 . 0 


ARE 


3 .3! 


ARE' 


ARE  +1.0 


2.0  -  6.0 


ARE' 
4 


ARE 


( AREq  +  1.0)' 


+2.0  -  (3.667  *  ARW3  +  2  *  ARW  +  1.667)*  Kp 

.2 

(2.5.18) 


20.0  *  ARE' 
(ARE3  +  2 )2 


TAN( LESWEEP)  represents  tne  approximate  value  of  BETA  at  which  CDW 
will  maximize,  provided  the  body  is  essentially  cylindrical  where 
the  wing  is  attached.  If  the  body  is  area-ruled,  the  peak  value 
of  CDW  may  or  may  not  be  closely  approximated. 

2 . 5 . 4 . 2  Body  Wave  Drag 

The  fuselage  body  wave  drag  is  computed  by  dividing  the  body 
into  two  parts,  consisting  of  a  simplified  pointed  nose  and  a 
simplified  boattail.  That  is, 

CDW  =  ( CDPN  +  CDBT )  *  AMAX/SREF  (2.5.19) 


Nose  wave  drag,  CDPN,  is  determined  from  Linnell's  empirical 
equation  for  the  supersonic  wave  drag  of  parabolic  noses  (Ref. 

18). 

CDPN  =  (1.2  +  1.5  *  X ) / (  (  1  . 0  +  1.9  *  X)  *  (1.0  +  ELOD2  )  ) 

(2.5.20) 


The  nose  fineness  ratio  ELOD  is  calculated  from  the  nose  length, 
LNOS  and  the  maximum  cross-sectional  area,  AMAX,as 

ELOD  =  LNOS/SQRT( AMAX  *  4.0/PI)  (2.5.21) 
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Boattail  wave  drag,  CDBT  is  determined  as  a  function  of  the 
boattail  fineness  ratio,  ELOD,  base  diameter  to  maximum  diameter, 
DBOD,  and  Mach  number.  This  is  done  by  computing  CDBT  at  five 
values  of  DBOD  and  interpolating  to  the  desired  value.  The 
general  form  of  these  equations  is  given  below: 

For  X  £  1 

CDBT(i)  =  ( ( i )  +  A,(i)  *  X  +  A,(i)  *  X2  +  A, ( i )  *  X3)/ELOD2 

°  (2.5.22) 


and  for  X  >  1 

CDBT(i)  =  A4 ( i ) / ( X  *  ELOD2)  (2.5.23 

where 

X  =  BETA/ELOD 

The  polynomial  coefficients  of  equations  2.5.22  and  2.5.23 
are  determined  from  a  least-square  fit  of  Fig.  III.B.10-9  of 
Aerospace  Handbook  (17)  for  ogive  boattails  and  are  tabulated 
below: 


I 

DBOD 

A0 

A1 

A2 

_ ^3 _ 

_ A4_ 

1 

0.0 

1.165 

-0.5112 

-0.5376 

0.3964 

0.513 

2 

0.4 

1  .  067 

-1  .709 

1 .6632 

-0.686 

0.3352 

3 

0.6 

0.7346 

-1.4618 

1 . 5795 

-0.6542 

0.198 

4 

0.8 

0.2555 

-0.5008 

0 .5024 

-0.2077 

0.0494 

5 

1  .  0 

0.0 

0.0 

0.0 

0.0 

0.0 

2.5.4. 3  Nacelle  Wave  Drac 


The  nacelle  wave  drag  is  calculated  by  a  method  similar  to 
that  used  for  the  fuselage: 


CDW  =  ( CDON  +  CDBT)  *  AMAX/SREF 


(2.5.24) 


where  CDBT,  AMAX,  and  SREF  are  as  before  and  the  equation  used  to 
calculate  CDON  for  open-nose  bodies  is 

CDON  =  ((1.0  -  2.0  *  RIN/DMAX ) /ELOD ) 1  ' 3/SQRT( BETA  )  (2.5.25) 

where 

RIN  is  the  nose  inlet  area  and  DMAX,  ELOD,  and  BETA  are  as  before. 

This  equation  is  a  curve  fit  of  Figure  III.B.10-6  of 
Aerospace  Handbook  (17). 


I 


2.5.5  Camber  Drag 

The  minimum  drag  contribution  of  the  surface  twist  and  camber 
is  related  to  the  lift  coefficient  of  the  polar  displacement, 
POLARD,  by  the  equation 


CDC  =  1  -  0/ ( PI 


(1.0  -  E) ) 


POLARD 


(2.5.26) 


This  increment  is  called  camber  drag  and  represents  a  drag 
increment  between  minimum  profile  drag  and  drag  minimum.  The  span 
efficiency  value,  E,  is  related  to  the  induced  drag  factor  and  is 
computed  by 


E  =  SQRT( EAR  *  ELE ) 

EAR  =  1.0  -  .45  *  ARS ‘ 68 

ELE  =  1  .  0  -  .000392  *  ABS ( LESWEEP ) 1  *  6 1 5 


(2.5.27) 


If  for  some  reason  E  >  1,  and  alternate  equation,  obtained  from 
The  Aerospace  Handbook  (17),  is  used 


CDC  =  .7  *  POLARD2  *  AEXP/ SERF 


(2.5.28) 


where 

AEXP 

ARS 

LEWSWEEP 


is  the  exposed  area  of  the  surface 

is  the  surface  aspect  ratio 

is  the  surface  leading  edge  sweep  angle 


2.5.6  Base  Drag 

Data  represented  in  Hoerner  (19)  were  used  to  establish 
equations  from  which  the  base  drag  of  bodies  could  be  determined. 
The  trends  of  these  data  show  three  different  phases:  (1)  a 
gradual  rise  in  CDB  at  transonic  speeds  up  to  Mach  =  1.0,  (2)  a 
relatively  constant  drag  level  supersonically  up  to  about  Mach  = 
1.8  and  (3)  a  steadily  decreasing  value  of  drag  above  Mach  =  1.8. 
The  resulting  empirical  equations  are  given  as 

CDB  =  (.1  +  .122  *  MACH8)  *  SBASE/SREF;  WHEN  MACH  £  1 

CDB  =  .222  *  SBASE/SREF;  WHEN  1  <  MACH  £  1.8 

CDB  =  (1.42  *  SBASE/SREF )/( 3. 1 5  +  MACH2);  WHEN  MACH  >  1.8 

(2.5.29) 

where 

SBASE  is  the  base  area. 


2.5.7  Drag  Rise 

For  Mach  numbers  less  than  Mach  critical  the  drag  increases 
slowly  with  increasing  Mach  number.  This  drag  component  is  known 
as  compressible  drag,  or  drag  creep.  Methodology  for  estimating 
this  component  of  drag  for  conventional  or  supercritical  wings 
was  included  in  the  subsonic  drag  buildup  described  above.  For 
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Mach  numbers  greater  than  Mach  critical,  drag  rise  begins  and 
increases  rapidly  with  Mach.  Figure  2.5.10  illustrates  the  drag 
bookkeeping  system  followed  in  SACP.  This  is  the  same  as  the 
system  used  in  the  Large  Aircraft  program  whereby  beyond  Mach  1.0 
the  drag  rise  and  the  interference  plus  form  drag  are  replaced  by 
wave  drag.  The  drag  rise  is  separated  into  two  components,  drag 
rise  due  to  lifting  components  and  drag  rise  due  to  all  other 
components  on  the  aircraft.  However,  at  this  point  only  the  total 
drag  rise,  at  zero  lift  due  to  lifting  and  non-lifting  components 
is  addressed.  This  value  is  determined  by 

CDDR  =  A2  *(MACH  -  MCR ) 2  +  A3  *(MACH  -  MCR ) 3  (2.5.30) 

where  A2  and  A3  are  defined  to  produce  a  continuous  zero-lift  drag 
curve  between  MCR  and  MACH  1.0.  The  drag  rise  is  curve  fitted  to 
begin  at  MCR  with  zero  slope  and  end  at  Mach  =  1.0,  matching  the 
value  and  slope  of  the  wave  drag  curve.  The  coefficients  A2  and 
A3  are  calculated  from 

A2  =  (3.0  * ( CDW1  -  CDFF )  -  (1.0  -  MCR)  *  CDWP ) / ( 1 . 0  -  MCR)2 

A3  =  ((1.0  -  MCR)  *  CDWP  -  2.0  *  ( CDW1  -  CDFF))/ (1.0  -  MCR)3 

where  CDW1 ,  CDFF,  and  CDWP  represent  the  Mach  1.0  value  of  the 
wave  drag,  form  plus  interference  drag,  and  the  slope  of  the  wave 
drag  respectively. 


2.6  SURFACE  INTERACTIONS 

Interaction  between  various  aircraft  surfaces  has  an  impact 
on  several  computed  properties.  The  interactions  between  canards 
and  wings  or  between  wings  and  horizontal  tails  impact  the 
calculations  of  angle  of  attack  for  these  components,  and  also  can 
affect  the  lift  curve  slope  of  the  wing.  The  methods  described 
below  make  the  necessary  adjustments  to  these  values  to  account 
for  the  interactions. 

2.6.1  Downwash  Effects  on  Angle  of  Attack 

Air  passing  over  a  surface  causes  a  disturbance  in  the 
airflow,  much  like  the  wake  produced  by  a  boat.  This  disturbance 
effects  other  surfaces  that  fall  within  the  wake  area.  This 
phenomenon  is  known  as  Downwash.  Downwash  alters  the  angle  of 
attack  on  all  of  the  aircraft  surfaces.  A  set  of  equations  for 
the  calculation  of  the  downwash  has  been  adapted  from  the  AEROX 
Computer  Program  for  Transonic  Aircraft  (20  and  21)  for  use  in 
SACP.  These  equations  are: 

T 1  =  2.0/ARS 

T2  =  SQRT( 2 . 0  *  SPANS/ (3.0  *  MACS)) 

T3  =  SQRT(SPANL) 

T4  =  (1.3  -  CONAR ) / ( SQRT (1.0  +  TRS  )  ) 

SPANL  =  .5  *  SPANS/LTS  . 

CONAR  =  1.0/ (3.0  *  ARSARS  3' 

OSPNL  =  1.0/ SPANL 

The  wake  angle,  EPSD,  is 

EPSD  =  ALCAM  *  T4  *  (T1  +  .1  *(T2  +  T3 ) )  (2.6.1) 

ALCAM  =  ALPHA  +  CLO/ ( 2 . 0  *  PI  *  C4SWEEP  *  ARS/(ARS  +  2.0)) 

the  nondimensional  offset  distance  from  the  affected  surface  to 
the  wake,  AYOSS,  is 

AYOSS  =  ABS (((2.0  *  ZS) /SPANS)  +  YDOSS  )  (2.6.2) 

YDOSS  =  ((OSPNL  -  1.5/ARS)  *  SIN(EPSD))  -  (OSPNL  *  SIN( ALPHA)) 

W1  =  . 2/ ( SQRT( OSPLN  )  ) 

W2  =  (2.0  +  COS (2.25  *  ALR))/3.0 
W3  =  1.0  -  (1.5  *  AYOSS ) 

W4  =  W2  -  .75  *  ALPHA  *  SIN(2.25  *  ALPHA) 

the  downwash  angle  at  the  effected  surface,  DWASH,  is 


DWASH 


ALCAM  *  T4  *  (W1  +  T1 )  *  W2  W3  *  ATENF 


(2.6.3) 
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and  the  downwash  derivative  with  respect  to  angle  of  attack. 


DEPDA, 

is 

DEPDA  = 

T4  *  (W1  +  1 

where 

ATENF 

is 

the 

is 

only 

ATENF  = 

1.0- 

.  1  * 

and 

(T1  +  i 

OSPNL 

ARS 

is 

the 

SPANS 

is 

the 

MACS 

is 

the 

LTS 

is 

the 

of 
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TRS 

is 

the 

ALPHA 

is 

the 

ZS 

is 

the 

C4SWEEP 

is 

the 

LESWEEP 

is 

the 

T1 )  *  W4  *  W3  *  ATENF 


(2.6.4  ) 


(2.6.5) 


aspect  ratio  of  the  effecting  surface 
span  of  the  effecting  surface 

mean  aerodynamic  chord  of  the  effecting  surface 
X  distance  between  the  aerodynamic  centers 
effecting  and  affected  surface 
taper  ratio  of  the  affecting  surface 
total  aircraft  angle  of  attack 
Z  distance  between  the  two  surfaces 
quarter  chord  sweep  of  the  effecting  surface 
leading  edge  sweep  of  the  effecting  surface 


This  algorithm  was  developed  for  the  downwash  angle  ( DWASH ) 
downstream  of  the  wing  at  aft  tail  locations.  For  other 
interactions,  the  downwash  angles  must  be  adjusted.  For  use 
downstream  of  canards,  the  angle  must  be  multiplied  by  the  ratio 
of  the  canard  span  to  the  wing  span 


ECW  =  DWASH  *  SPANC/SPANW 


(2.6.6) 


For  canards  mounted  outboard  of  widely  spaced,  wing  mounted 
nacelles  (NCAN  =  2)  where  the  wing  inboard  panels  are  out  of  the 
canard  downwash, 


ECW  =  0 


The  wing  will  also  produce  an  upwash  disturbance  at  the 
canard.  In  this  case,  the  ratio  of  the  upwash  velocity  (upstream 
of  the  lateral  or  bound  leg  of  the  equivalent  "horseshoe"  vortex 
representing  the  wing)  to  the  downwash  velocity  at  an  equal 
distance  downstream  of  the  lateral  vortex  leg  was  evaluated  by  the 
Biot-Savart  law  for  vortex-induced  velocities.  This  ratio  is 
denoted  by  FWAC. 


FWAC  =  1.0  -  (1.2  *  (LTO/SPANW) 


(2.6.7) 
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The  upwash  angle  ahead  of  the  wing  is 
EWC  =  DWASH  *  FWAC 

Three  additional  conditions  are  now  placed  on  the  values  of  wing 
upwash  angle.  First,  for  sharp  wings  having  aspect  ratios  greater 
than  2.5,  the  leading-edge  Mach  limited  flow  zone  is  assumed  to 
cancel  the  upwash 

EWC  =  0;  WHEN  NATW  =  9  ( B 1  convex  airfoil  type)  and  ARW  >  2.5 

Second,  for  canards  mounted  outboard  of  widely  spaced,  twin 
nacelles  ( NCAN  =  2),  the  wing  upwash  angle  is  reduced  by  a  factor 
involving  the  spans  of  the  wing  and  canard  and  the  nacelle  spacing 
( SPANAC  )  . 

EWC  =  DWASH  *  FWAC  *  (1.0  -  (SPANAC  +  SPANC/ 2 . 0 ) / SPANW )  (2.6.8) 

Third,  for  supersonic  Mach  numbers,  the  wing  upwash  is  attenuated 
by  dividing  by  the  Mach  number. 


EWC  =  DWASH  *  FWAC/MACH ;  WHEN  MACH  >  1.0 


(2.6.9) 


The  relationship  of  the  wing  to  the  horizontal  tail  is  similar 
with : 

EWH  =  DWASH 

FHAW  =  1.0  -(1.2  * ( LTW  )  / SPANH ) 

EHW  =  DWASH  *  FHAW 

EHW  =  DWASH  *  FHAW/MACH ;  WHEN  MACH  >  1.0 

EHW  =  0;  FOR  SHAPE  LEADING  EDGE  AIRFOILS  AND  ARH.GT.2.5 

The  downwash  and  upwash  angles  produced  by  the  canard,  wing,  and 
horizontal  tail  can  now  be  used  to  obtain  the  angles  of  attack  for 
these  three  surfaces  including  the  induced  angle  of  attack: 


ALPHAC  =  ALPHA  +  IC  +  EWC 


(2.6.10) 


ALPHAW  =  ALPHA  +  IW  -  ECW  +  EHW  +  EWC  *  SPANC/2.0  *  SPANW) 

(2.6.11  ) 


ALPHAH  =  ALPHA  +  IH  -  EWH  +  EHW  *  (SPANW/2.0  *  SPANH) 


(2.6.12) 


where  IC,  IW,  and  IH  are  the  incidence  angles  of  the  canard,  wing, 
and  horizontal  tail. 

2.6.2  Wing  Lift 

The  possibility  of  favorable  interference  from  close-coupled 
canards  on  wing  lift  was  cited  initially  by  Behrhohm  (22).  An 
aerodynamic  lift  interference  factor  (VORTEX)  has  been  added  to 
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SACP  and  may  increase  the  upper  surface  lift  of  wings  having 
aspect  ratios  near  2.0  and  having  close-coupled  canards  mounted 
above  the  wing  plant.  This  factor,  in  effect,  accounts  for  the 
"fluid  end-plate  action"  of  the  canard  trailing  vortices  which 
promote  local  two-dimensional  flow  over  the  inboard  wing  upper 
surface,  when  the  wing  operates  with  a  subsonic  leading  edge.  No 
evidence  of  similar  vortex  action  has  been  observed  with  wings 
having  supersonic  leading  edges.  The  interference  factor  is 
proportional  to  the  square  root  of  the  absolute  value  of  the 
canard  lift  coefficient  (ABCLC),  and  decreases  with  further 
upstream  positioning  of  the  canard  ( GLC )  and  with  vertical 
repositioning  of  the  canard  above  or  below  the  favorable,  close 
coupled  height  zone  (GCZ). 

GCL  =  1.4  *  MACW/LTC ;  WHEN  LTC/MACW  _>  1  . 4 

GCL  =  1.0;  WHEN  LTC/MACW  <  1.4 

GCZ  =  1.0;  WHEN  0.0  <  ZCAN/MACW 

GCZ  =  1.0  +  ZCAN/MACW;  WHEN  0.0  >  ZCAN/MACW 

GCZ  =  1.0  -  ZCAN / ( 4 . 0  *  MACW);WHEN  .2  <  ZCAN/MACW 

For  1.75  <  ARW  <  2.5,  the  interference  factor  is 

VORTEX  =  GCL  *  GCZ  *  CLC/(ARWC  *  ABCLC) 

*  SQRT( ABCLC  *  SPANC  *  MACW/ ( SPANW  *  MACC )  ) 

*  (1.0  -  ABS ( ARW  -  2  .  1  25  )/.  375)  (  2.6.1  3) 

where 

ARWC  =  ARW/ (ARW  +  2.0) 

CLC  is  the  canard  lift  coefficient  referenced  to  the  wing 

area 

MACW  is  the  mean  aerodynamic  chord  of  the  wing 

MACC  is  the  mean  aerodynamic  chord  of  the  canard 

LTC  is  the  X  distance  between  canard  and  wing  aerodynamic 

center 

ZCAN  is  the  Z  distance  between  canard  and  wing  aerodynamic 
center 

ARW  is  the  aspect  ratio  of  the  wing 

SPANC  is  the  span  of  the  canard 

SPANW  is  the  span  of  the  wing 


The  interference  factor  is  reduced  for  wings  having  highly 
swept  leading  edges,  because  some  leading  edge  vortex  lifts 
already  act  on  the  wing  with  or  without  canard.  For  wings  having 
leading  edge  sweep  angle  greater  than  45°  and  aspect  ratio  between 
1.75  and  2.5, 

VORTEX  =  VORTEX/TAN ( LESWEEP  ) 


For  supersonic  Mach  numbers  and  subsonic  wing  leading  edges, 
VORTEX  =  VORTEX/MACH 

The  lift  coefficient  of  the  wing  upper  surface  and  the  wing  lift 
curve  slope  can  now  be  computed  by  multiplying  the  interference 
factor  as  follows: 

CLWU  =  CLWU  *(1.0  +  VORTEX )  (2.6.14) 

CLWALPHA  =  CLWALPHA  *  (1.0  +  VORTEX/2.0)  (2.6.15) 

Again,  the  interference  factor  is  used  only  for  wings  having 
aspect  ratios  between  1.75  and  2.50. 


2.7  LIFT,  LIFT  CURVE  SLOPE,  AND  DRAG  DUE  TO  LIFT 
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Methods  for  the  computation  of  lift,  lift  curve  slope  and 
drag  due  to  lift  were  taken  from  the  AEROX-Computer  Program  (20 
and  21)  for  use  in  SACP.  The  equations  are  organized  according  to 
the  division  of  the  flight  envelope  into  the  flow  zones  depicted 
in  figure  2.7.1.  The  cross-hatched  viscous  stall  region  is  not 
included  in  the  present  method.  The  flow  is  considered 
incompressible  for  Mach  numbers  below  0.1.  The  compressible, 
shockless  zone  2  covers  airfoils  with  blunt  leading  edges  and  with 
surrounding  flows  having  local  Mach  numbers  everywhere  below  the 
Laitone  limit  value  (SQRT(1.4  +  3.0)/2.0).  The  boundary  between 
zones  2  and  3  is  estimated  by  the  equation  for  ALER,  the  angle  of 
attack  for  onset  of  the  limit  Mach  number  at  the  leading  edge. 

The  onset  of  the  surface  limit-Mach  number,  zone  4,  occurs  when 
the  supper-surface  lift  for  zone  2  or  3  reaches  the  limit  lift 
corresponding  to  the  designated  chordwise  shock  location,  XCD. 
Sharp  airfoils  are  considered  to  have  reached  the  leading-  edge 
Mach  limit  for  all  subsonic  flow,  and  are  treated  with  the 
nonpotential  lift  equations  of  zones  3  and  4.  Zone  5  applies  only 
to  sharp  airfoils  having  supersonic  leading  edges  with  attached 
shocks.  Finally,  Flow  zone  6  covers  all  airfoils  with  detached, 
leading  edge  shocks.  The  supersonic  leading  edges  for  zones  5  and 
6  occur  when  the  normal  component  of  Mach  number  exceeds  unity 
(Mach  *  COS  ( LESWEEP )  1.0). 

The  following  equations  apply  to  an  equivalent  wing  having 
straight  leading  and  trailing  edges.  Account  is  made  for  strakes 
and  wing  forward-chord  extensions  in  SACP  by  multiplying  all  of 
the  following  equations  for  lift  coefficient  and  wing  lift  curve 
slope  by  the  term  (1.0  +  FLEX).  FLEX,  the  empirical  lift  factor 
for  chord  extension,  is  defined  within  the  program  and  depends  on 
Mach  number,  angle  of  attack,  and  whether  the  extension  is  sharp 
( IFLEX  =  0 )  or  blunt  ( IFLEX  =1). 

2.7.1  Incompressible  Flow  Zone  1  (M  <  0.1) 


The  incompressible  lift  equation  for  conventional  airfoils  is 
taken  from  the  potential-flow  theory  of  Kutta-Jowkowski .  The  lift 
equation  for  nonpotential  flow,  such  as  around  sharp  airfoils,  is 
based  upon  the  integration  of  downwash  momentum.  Both  equations 
are  extended  to  three-dimensional  flows  through  the  inclusions  of 
the  Prandtl  aspect-ratio  transformation  and  the  cosine  term 
involving  the  effective  sweep  of  the  quarter-chord  line.  The 
potential-flow  equations  for  conventional  airfoils  are: 

CL  =  2.0  *  PI  *  SIN( ALPHA)  *  COS(SQ)  *  AR/(AR  +2.0)  (2.7.1) 

DLL  =  2.0  *  PI  *  C0S( ALPHA)  *  COS(SQ)  *  AR/(AR  +2.0)  (2.7.2) 

CD  =  CL2/( PI  *  AR)  (2.7.3) 


TOW 


MACH  NUMBER 


where 
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ALPHA  is  the  surface  angle  of  attack 
AR  is  the  surface  exposed  aspect  ratio 

SQ  is  the  effective  quarter  chord  sweep  of  the  surface 

SQ  =  ARLSIN  ( SIN ( SWEEPC4 )  *  COS (ALPHA)) 

The  nonpotential  equations  derived  for  AEROX  (20)  are: 

CL  =  2  *  PI  *  SIN( ALPHA )  *  COS(ALPHA)2  *  (1  -  SIN( ALPHA ) 2^2 )  * 

COS(SQ)  *  AR/ ( AR  +2)  (2.7.4) 

DCL  =  PI  *  C0S( ALPHA)  *  (2  *  SIN( ALPHA)4  +  2  *  C0S( ALPHA)2  - 

4  *  SIN ( ALPHA ) 2  -  3  *  SIN ( ALPHA ) 2  *  COS(ALPHA)2)  *  COS(SQ) 

*  AR/ ( AR  +2)  (2.7.5) 

=  PI  *  COS( ALPHA)  *  (CISQ)  *  COS(SQ)  *  AR / ( AR  +  2) 

CD  =  CL  *  TAN ( ALPHA )DEXP  (2.7.6) 

where 

DEXP  =  1.5  -  MACH  *  COS ( LESWEEP ) / 4 ;  WHEN  M  <  1 
=  1.5  -  MACH2  *  COS ( LESWEEP )/4;  WHEN  M  >  1 

The  incompressible  lift  is  comprised  of  equal  contributions 
from  the  upper  and  lower  wing  surfaces.  The  nonpotential 
equations  asymptotically  approach  the  potential  equations  at  small 
angles  of  attack. 

2.7.2  Compressible,  Shock-free  Flow  Zone  2 

In  estimating  the  aerodynamics  for  airfoils  having  subsonic 
leading  edges,  separate  compressibility  factors  are  used  for  the 
lifts  of  the  upper  and  lower  surfaces.  For  compressible,  shock- 
free  flow  (zone  2)  on  the  upper  surface  of  blunt  airfoils,  a 
Prandt 1-Glauert  factor  is  used  involving  the  component  of  flight 
Mach  number  normal  to  the  quarter-chord  line  and  an  angle-of- 
attack  attenuation  to  value  unity  at  40°. 

FU  =  1.0  -  (1.0  -  SQRTd.O  -  MACH2  *  COS(SQ)2)  *  (ALPHA/40)2/ 

SQRTM.O  -  MACH2  *  COS  (  SQ  )  2  )  (2.7.7) 

No  compressibility  factor  is  used  for  the  upper-surface  lift 
in  nonpotential  flow,  because  the  onset  of  the  local  Mach-limit  is 
considered  to  "freeze"  the  local  flow. 

Lower  surface  lifts  in  potential  and  nonpotential  flows 
around  airfoils  having  subsonic  leading  edges  are  evaluated  using 
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a  pitot-compressibility  factor  expressing  the  ratio  of  the 
leading-edge  stagnation -line  pressure  in  compressible  flow  to  that 
in  incompressible  flow.  Account  is  included  for  variations  in  the 
angles  of  attack  and  sweep,  and  for  wing-body  interference 
patterned  after  the  FLAX  parameter  from  references  2  4  and  25. 

FL  =  CPMXS  =  1  . 4^8  57  1  7/MACH2  *  ((  1.0  +  .2  *  FINT  *  MACH2  * 

C0S( ESWPLE)^  )3 • 5  -  1  . 0 VC0S( ESWPLE )2  (2.7.8) 

where 

ESWPLE  is  the  effective  leading  edge  sweep  of  the  surface 

=  ARCSIN  ( SIN( LESWEEP  )  *  C0S( ALPHA)) 

FINT  =  ((1.0  +  BDMAX ) /SPANS )2 ;  WHEN  1.4  <  MACH 

FINT  =  (1.0  +  (1.0  -  (ALPHA/40)  *  (1.4  -  MACHJ/.04)2  * 

BDMAX/ SPANS  ) 2 ;  WHEN  (1  <  MACH  <_  1.4) 

FINT  =  (1.0  +  (1.0  -  ALPHA/40.)2  *  BDMAX/ SPANS ) 2 ;  WHEN  (MACH  £  1) 

BDMAX  is  the  body  width 

SPANS  is  the  span  of  the  surface 

The  compressible,  potent ia 1 - f low  equations  for  conventional 
airfoils  are  obtained  by  incorporating  the  compressibility 
factors,  equations  2.7.7  and  2.7.8  into  equations  2.7.1  and  2.7.2. 

CL  =  (FU  *  FINT  +  FL)  *  PI  *  SIN (ALPHA)  *  COS(SQ)  *  AR/(AR  +2.0) 

(2.7.9) 

DCL  =  ((FU  *  COS ( ALPHA )  +  ( dFU / dALPH A )  *  SIN(ALPHA))  *  FINT  + 

FL  *  C0S( ALPHA))  *  PI  *  COS(SQ)  *  AR/(AR  +  2.0)  (2.7.10) 

CD  =  CL2/ (PI  *  AR) 

Sharp  airfoils  in  subsonic,  compressible  flow  are  treated  as 
nonpotential,  Mach-limited  flows,  and  assigned  to  flow  zones  3  or 

4  . 

2.7.3  Leading-edge  Mach-lim.ited  Flow  Zone  3 

The  onset  of  the  Laitone  limit  Mach  number  in  the  curvilinear 
flow  around  a  blunt-air foi 1  leading  edge  is  estimated  by  the  value 
of  the  onset  angle  of  attack,  ALPHAE.  A  brief  discussion  of  the 
conceptual  flow  model  undergoing  the  transition  from  compressible, 
potential  flow  to  Mach - limited ,  ncnpotent i a  1  flow  follows.  For 
more  information  see  the  appendix  in  AERC'X  (20). 


After  attaining  the  local  limit  Mach  number  around  the 
leading  edge,  further  increases  in  angle  of  attack  or  in  flight 
Mach  number  tend  to  enlarge  the  local  Mach-limited  region  on  the 
airfoil  nose. 

The  enlargement  of  the  constant  velocity  profile  in  the 
curvilinear  flow  around  the  nose  disrupts  the  radial  equilibrium 
between  the  local,  radial  static-pressure  gradient  and  the  local, 
centrifugal  focus.  ( Irrotational ,  curvilinear  flow  exists  only 
when  radial  equilibrium  prevails,  which  requires  an  essentially 
inverse  relationship  between  local  velocity  and  local  streamline 
radius  of  curvature.)  The  constant  velocity  profile  may  be 
accomplished  by  the  dominance  of  the  local  centrifugal  forces  and 
the  possible  formation  of  a  separation  bubble.  When  the  locally 
supersonic  flow  decelerates  (and  reattaches)  just  downstream  of 
the  nose  (with  or  without  the  separation  bubble),  a  "peaky" 
pressure  distribution  results,  (shown  in  ref.  26).  This  flow  is 
designated  in  AEROX  as  leading-edge  Mach-limited  or  zone  3.  When 
the  separation  bubble  is  present,  the  mixing  action  of  the 
rotational-flow  lagers  emanating  from  the  disrupted  nose  flow  may 
promote  the  flow  reattachment. 

When  the  supersonic  limit-Mach  number  extends  well  back  on 
the  airfoil  upper  surface,  resulting  in  the  flat  pressure 
distribution  (also  in  ref  26),  the  flow  is  surface  Mach-limited, 
and  assigned  to  zone  4.  Surface  Mach-limited  flow  is  modeled  to 
have  separation  downstream  of  the  surface  limit  shock,  resulting 
in  an  additive  drag  component,  CDSEP.  Experimental  confirmation 
of  the  separation  appears  in  pressure  distributions,  samples  of 
which  were  included  in  reference  26. 

Sharp  airfoils  in  compressible,  subsonic  flow  are  treated  as 
leading-edge  Mach-limited  flow  with  the  nonpotential  zone  3 
equations . 

CL  =  (1.0  +  FL)  *  PI  *  SIN ( ALPHA )  *  COS(ALPHA)2  *  (1  - 

SIN ( ALPHA ) 2 / 2 . 0 )  *  COS(SQ)  *  AR/(AR  +  2.0)  (2.7.11) 

DCL  =  (1  +  FL)  *  PI  *  C0S( ALPHA)  *  (CISQ)  *  COS(SQ)  *  AR/(AR  +  2) 

(2.7.12) 

CD  =  CL  *  TAN( ALPHA )DEXP  (2.7.13) 

The  lift  curves  for  blunt  airfoils  traversing  from 
compressible,  potential  flow  (zone  2)  into  zone  3  are  assigned  a 
continuous  mathematical  transition,  rather  than  a  discontinuous 
jump,  so  that  the  program  may  be  coupled  to  an  optimizer  program. 
The  lift  equation  is  the  sum  of  the  zone  2  lift  (eg.  2.7.9) 
evaluated  at  ALPHAE  and  an  incremental  lift  for  the  interval 
(ALPHA  -  ALPHAE)  in  zone  3,  until  equation  2.7.11  is  reached. 
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CL  =  The  greater  of 
Equation  2.7.11 


CL  =  (FU  *  FINT  +  FL)  *  PI  *  SIN ( ALPHAE )  *  COS(SQ)  *  AR/(AR  + 

2.0)  +  (dCLU/dALPHA  +  dCLl/dALPHA )( ALPHA  -  ALPHAE)  (2.7.14) 

2.7.4  Surface  Mach-Limited  Zone  4 

The  conceptual  flow  model  depicting  surface  Mach-number 
limited  flow  past  an  airfoil  is  shown  in  figure  2.7.2.  The 
Laitone  limit  Mach  number  (ref.  27  and  28)  corresponds  to  the 
maximization  of  the  static  pressure  behind  the  surface  limit 
shock . 


MACHlim  =1.48 


(2.7.15 


This  criterion  is  extended  to  swept  wings  in  the  AEROX 
program  with  the  following  limit  pressure  coefficient. 


CPLIM  =  1  .4285717/MACH2  *  ( 1.0  +  .2  *  MACH2  *  COS ( SXC  ) 2  ) 3 * 5/ 


where 


3.58  -  1.0 


is  the  effective  sweep  angle  of  the  limit  shock 


(2.7.16 


The  airfoil  lift  coefficient  for  zone  4  is  the  sum  of  the 
lower-surface  nonpotential  lift  (equation  2.7.11)  and  the  upper- 
surface  limit  lift,  which  is  dependent  upon  the  designated 
chordwise  location,  XCD  of  the  surface  limit  shock.  The  program 
uses  XC  =  f(XCD)  which  moves  the  limit  shock  to  the  trailing  edge 
as  the  sonic  leading-edge  condition  is  reached. 

CL  =  -(CPLIM  *  XC  +  .5  *  CPDS  *  (1.0  -  XC ) )  *  C0S( ALPHA)  + 

FL  *  PI  *  SIN(ALPHA)  *  COS ( ALPHA ) 2  *  (1  -  SIN ( ALPHA  ) 2 / 2 )  * 
COS(SQ)  *  AR/(AR  +2)  (2.7.17) 

DCL  =  (CPLIM  *  XC  +  . 5 ( CPDS  )  *  (1.0  -  XC )  )  *  SIN( ALPHA )  +  FL  * 

PI  *  C0S( ALPHA)  *  ( CISQ )  *  CCS(SQ)  *  AR/(AR  +  2.0)  (2.7.18) 


CD  =  CL2/ (PI  *  AR )  +  CL  *  TAN( ALPHA )DEXP  +  CDSEP 


(2.7.19) 


where 


CPDS  is  the  pressure  coefficient  downstream  of  the  limit  shock, 
corresponding  to  P2  in  figure  2.7.2 
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CDSEP  is  the  separation  drag.  It  accounts  for  the  momentum 

deficit  in  the  modeled  separation  wake  passing  the  trailing 
edge.  The  wake  is  assumed  to  originate  at  the  base  of  the 
limit  shock,  to  have  its  upper  edge  follow  a  line  inclined 
at  one-half  the  angle  of  attack,  and  to  have  a  linear 
velocity  profile  between  the  zero  value  on  the  surface  to 
the  freestream  value  at  the  upper  edge. 

CPDS  =  1.4285714/MACH2  *  (((2.8  *  (2.2  +  MACH2  *  SIN(SXC)2)  * 
COS(SXC)2  -  .  4  )  / 2 . 4  )  *  ((2.  +  .4  *  MACH2)/  (.4  *  (2.2  + 

MACH2  *  SIN ( SXC  )  2  +  2 .  )  )  )  3  *  5  -  1.) 


CDSEP  =  (SPANS  -  BDMAX )  *  (l.  -  XC 5  *  MACS/SREF  *  SIN(ALPHA)/2 


where 

SPANS 

BDMAX 

MACS 

SREF 

SXC 


is  the  span  of  the  surface 
is  the  body  diameter 

is  the  mean  aerodynamic  chord  of  the  exposed  surface 
is  the  reference  area 

is  the  effective  sweep  angle  of  the  limit  shock 


2.7.5  Supersonic  Attached-Shock  Zone  5 


Airfoils  having  sharp  leading  edges  with  attached  shocks  are 
treated  in  zone  5.  The  upper  boundary  is  the  angle  of  attack  for 
shock  detachment,  ARDET,  which  is  evaluated  by  equations  based  on 
curve  fitting  figure  4  in  reference  29  for  supersonic  speeds.  The 
supersonic  lift  coefficients  are  estimated  by  the  nonpotential 
equation  (2.7.11)  using  the  exposed  wing  area,  an  empirical,  Mach- 
attenuated  aspect-ratio  transformation,  and  the  wing-body 
interference  factor  of  FLAX  (ref.  24  and  25)  applied  to  the  lower 
surface  lift. 

CL  =  2.  *  PI  *  SIN( ALPHA )  *  COS(ALPHA)2  *  (1.  -  SIN ( ALPHA ) 2 / 2 .  * 
ARX/ ( ARX  +  2. /MACH)  *  SX/SREF  *  FLAX  *  (1.  -  .65/MACH)  + 
.65/MACH);  when  (MACH  <_  SQRT(2.)  (2.7.20) 

DCL  =  2 .  *  PI  *  COS ( ALPHA )  *  (CISQ)  *  ARX/(ARX  +  2 . /MACH  ) 

SX/SREF  *  FLAX ( 1 .  -  .65/MACH)  +  .65/MACH); 

when  (MACH  £  SQRT( 2 .  )  (2.7.21) 

For  SQRT(2.)  <  Mach  £  3.  multiply  the  above  equations  by 
1  . /SQRT( MACH2  -  1  .  ) 

At  MACH  numbers  above  3.,  lift  coefficients  are  evaluated  by  the 
explicit,  oblique-shock  theory  of  reference  30. 
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CL 


=  .833  *  SX/SREF  *  ARX/ ( ARX  +  2 . /MACH )  *  (1.  +  1.4  * 

SIN( ALPHA ) 2  -  COS ( ALPHA )  *  SQRT( 1 .  -  4 . /MACH2  *  ALPHA/16.  - 
1.96  *  SIN( ALPHA)2 ) )  *  COS( ALPHA); 

when  (ALPHA  16.  degrees;  and  3.  <  MACH)  (2. 7. 22. A) 

CL  =  .833  *  SX/SREF  *  ARX/ (ARX  +  2 . /MACH )  *  (1.  +1.4  * 

SIN ( ALPHA ) 2  -  COS ( ALPHA )  *  SQRT( 1 .  -  4. /MACH2  -  1.96  * 

SIN( ALPHA)2 ) )  *  COS (ALPHA);  when  (16.  degrees  £  ALPHA;  and 
3.  <  MACH)  (2.7.22.B) 

DCL  =  SX/SREF  *  (.833  *  2.8  *  SIN ( ALPHA )  *  COS(ALPHA)2  + 

SIN( ALPHA)  *  COS( ALPHA)  *  SQRT( 1 .  -  4. /MACH2  -  1.96  *  SIN 
(ALPHA)2)  +  (1.96  *  SIN( ALPHA)  *  COS ( ALPHA ) 3  ) / 

SQRT(1.  -  4. /MACH2  -  1.96  *  SIN( ALPHA ) 2 ) )  *  ARX/(ARX  + 

2. /MACH)  -  CL  *  TAN( ALPHA));  when  (3.  <  MACH)  (2.7.23) 

CD  =  CL  *  TAN( ALPHA)  (2.7.24) 

where 

ARX  is  the  aspect  ratio  of  the  exposed  surface 

SX  is  the  area  of  the  exposed  surface 

SREF  is  the  reference  area 

2.7.6  Supersonic  Detached-Shock  Zone  6 


Blunt  airfoils  with  supersonic  leading  edges  and  sharp 
airfoils  with  detached  leading  edge  shocks  are  treated  with  the 
following  equations  combining  the  nonpotential  equation  (2.7.11) 
derived  from  diverted  momentum  and  modified  Newtonian  impact 
theory,  which  assumes  increased  importance  on  the  windward  surface 
loading  at  hypersonic  speed.  The  limit  pressure  coefficient  for 
the  upper  surface  undergoes  a  transition  from  a  value  of  -.68  at 
MACH  =  1.  and  approaches  the  value  -l./MACH2  the  well-known 
empirical  value  suggested  by  Mayer  in  reference  31. 

CPLIM  =  -l./MACH2  *  (1.  -  .32/MACH2'5)  (2.7.25) 

The  upper  surface  lift  coefficient  becomes  a  small  fraction 
of  the  total  lift  coefficient  as  the  Mach  number  increases. 

CLU  =  the  lesser  of 

-CPLIM  *  SX/SREF  *  COS( ALPHA) 


or 

PI  *  SIN ( ALPHA )  *  COS ( ALPHA  )  2  *  (1.  -  SIN ( ALPHA  )  2 / 2 .  )  * 
SX/SREF  *  ARX/ (ARX  +  2 . /MACH  )  *  FMOMU/DNBL  (2.7.26) 


CLL  =  PI  *  SIN ( ALPHA )  *  COS(ALPHA)2  *  (1.  -  SIN( ALPHA ) 2 /2 . )  *  SX/ 
SREF  *  ARX/ ( ARX  +  2. /MACH)  *  FMOML/DNBL  *  FLAX  *  CPSTAG  + 
CPSTAG  *  SIN ( ALPHA ) 2  *  COS (ALPHA)  (2.7.27) 

CL  =  CLU  +  CLL 

DCLU  =  the  lesser  of 

CPLIM  *  SX/SREF  *  SIN( ALPHA) 

or 

PI  *  COS( ALPHA)  *  CISQ  *  ARX/ (ARX  +  2 . /MACH )  *  SX/ 

SREF  *  FMOMU/DNBL  (2.7.28) 

DCLL  =  PI  *  COS (ALPHA)  *  (CISQ)  *  ARX/ (ARX  +  2. /MACH)  *  SX/SREF  * 
FMOMU/DNBL  *  FLAX  *  CPSTAG  +  CPSTAG  *  SIN( ALPHA)  * 

(2.  *  COS ( ALPHA ) 2  -  SIN( ALPHA )2 )  (2.7.29) 

DCL  =  DCLU  +  DCLL 

CD  =  CL  *  TAN (ALPHA)  (2.7.30) 


The  empirical  constants  entering  the  diverted-momentum 
(nonpotential)  lift  contribution  include  DNBL ,  establishing  the 
variation  with  Mach  number,  and  FMOMU  and  FMOML,  expressing  the 
division  of  the  diverted  momentum  lift  to  the  upper  and  lower  wing 
surfaces . 

DNBL  =  MACH2/3;  when  MACH  <.1.77 

or  9 

SQRT( MACH^  -  1.);  when  MACH  >  1.77 

FMOMU  =  1  .  ;  when  MACH  <^1.2 

or 

1 . / ( MACH  -  2.);  when  MACH  >  1.2 

FMOML  =  1.6  -  .6  *  MACH;  when  MACH  (  1.5 

.7;  when  1 . 5  <  MACH  ^2.2 

or 

0.1  +  1 . / ( MACH  -  .55); when  MACH  >  2.2 


2.7.7  Lift  and  Drag  for  Bodies  and  Nacelles 


The  lift  and  drag  forces  on  the  body  or  on  the  nacelles  are 
estimated  by  the  equation  summing  the  slender  body  contribution 
and  the  viscous  cross  flow  drag.  An  updated  discussion  of  the 
approach  appears  in  reference  32,  which  was  the  principle  source 
of  the  values  used  for  constructing  the  contour  plot  of  cross- 


flow  drag  coefficient  against  cross-flow  Reynolds  number  and 
cross-flow  Mach  number  shown  in  figure  2.7.3.  The  BAERO 
subroutine  contains  explicit  equations  for  each  of  the  nine 
regions  indicated. 

CL  =  ( SIN( 2 .  *  ALPHA)  *  COS ( ALPHA/2 . )  *  ASECT  +  CDC  *  ETAN  * 

SIN ( ALPHA ) 2  *  APLAN ) /SREF  *  COS (ALPHA)  +  CLO  (2.7.30 


CD  =  (SIN(2.  *  ALPHA)  *  COS ( ALPHA/2 .  )  *  ASECT  +  CDC  *  ETAN  * 

SIN(ALPHA)2  *  APLAN) /SREF  *  SIN( ALPHA)  +  CDO  (2.7.31) 

where 

ETAN  =  .55  +  .0125  *  LENGTH/DIAMETER 

ASECT  is  the  nose  maximum  cross  sectional  area 

APLAN  is  the  nose  planform  area 

CLO  is  the  zero  AoA  lift  coefficient 

CDO  is  the  zero  AoA  drag  coefficient 


2.8  DYNAMIC  PRESSURE  RATIO 

The  dynamic  pressure  ratios  are  computed  for  horizontal 
downstream  surfaces.  In  addition  the  sidewash  times  the  dynamic 
pressure  of  the  vertical  tail  can  also  be  computed. 


2.8.1  Dynamic  Pressure  at  Down  Stream  Surfaces 


The  method  for  estimating  the  dynamic-pressure  QOQ  at  the 
downstream  surfaces  is  based  on  the  DATCOM  method  which  relates 
the  dynamic-pressure  ratio  to  the  drag  coefficient  of  the  upstream 
surface.  The  steps  involved  in  determining  the  dynamic  pressure 
at  some  distance  aft  of  the  wing  root  chord,  outlined  in  section 
4.4.1  of  the  DATCOM  report  are  as  follows: 

1 .  Calculate  the  half-width  of  the  wing  wake  by 

ZWOC  =  0.68  *  SQRT( CDOUS  *  ( XOC  +.15))  (2.8.1) 

where 

CDOUS  is  the  wing  zero-lift  drag  coefficient  of  the 

upstream  surface 

XOC  is  the  distance  measured  from  the  upstream  surface 

trailing  edge  to  the  down  stream  surface  leading 
edge 

2.  Calculate  the  downwash  in  the  plane  of  symmetry  at  the  vortex 
sheet  by 

EW  =  0.51566  *  (2.  *  CLALPKA )  / AR  (2.8.2) 

where 

CLALPHA  is  the  lift  curve  slope  of  the  upstream  surface 

AR  is  the  aspect  ratio  of  the  upstream  surface 

3.  Determine  the  vertical  distance  Z  from  the  vortex  sheet  to 
the  quarter-chord  point  of  the  mean  aerodynamic  chord  of  the 
downstream  surface  by 

ZOC  =  XOC  *  TAN( GAMMA  +EW-2.)  (2.8.3) 

where 

GAMMA  =  ATAN( ZMAC/XDIST ) 

ZMAC  is  the  Z  location  of  the  downstream  mean  aerodynamic 
chord 

XDIST  is  the  distance  in  feet  between  the  upstream  surface 
trailing  edge  and  the  downstream  surface  leading  edge 
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4.  Determine  the  dynamic-pressure-loss  ratio  at  the  wake  center 
by 

DQOQO  =  2.42  *  SQRT ( CDOUS ) / ( XOC  +  .3)  (2.8.4) 

5.  Determine  the  dynamic-pressure-loss  ratio  for  points  not  on 
the  wake  centerline  by 

DQOQ  =  DQOQO  *  (COS (PI/2  *  ZOZW))2  (2.8.5) 

where 

zozw  =  zoc/zwoc 

6.  Determine  the  dynamic  pressure  ratio  at  an  arbitrary  distance 
X  aft  of  the  upstream-surface-root-chord  trailing  edge  by 

QOQ  =  1  -  DQOQ 


2.8.2  Side-Wash-Factor  *  Dynamic  Pressure  Ratio  of  the  Vertical 
Tail 


For  a  wing-body  combination  there  are  two  contributions  to 
the  sidewash  present  at  a  vertical  tail  -  that  due  to  the  body  and 
that  due  to  the  wing. 

The  sidewash  due  to  a  body  arises  from  the  side  force 
developed  by  a  body  in  yaw.  As  a  result  of  this  side  force  a 
vortex  system  is  produced,  which  in  turn  induces  lateral-velocity 
components  at  the  vertical  tail.  This  sidewash  from  the  body 
causes  a  destabilizing  flow  in  the  airstream  beside  the  body. 

Above  and  below  the  fuselage,  however,  the  flow  is  stabilizing. 

The  sidewash  arising  from  a  wing  in  yaw  is  small  compared  to 
that  of  the  body.  The  flow  above  the  wake  center  line  moves 
inboard  and  the  flow  below  the  wake  center  line  moves  outboard. 


For  conventional  aircraft  the  combination  of  the  wing-body 
flow  fields  is  such  as  to  cause  almost  no  sidewash  effect  below 
the  wake  center  line. 


A  single  algebraic  equation,  taken  from  the  DATCOM,  is 
presented  that  predicts  the  combined  sidewash  and  dynamic  pressure 
parameter.  This  empirically  derived  expression  is 


SWFNV 


.724  +  (3.06  *((SV/SW)/1.  +  COS ( C4SWEEP ) ) )  + 
(.4  *  ZCGC4W/BMW ) )  +  (.009  *  ARW ) 


(2.8.6) 


where 


SV  is  the  vertical  tail  total  area 

SW  is  the  wing  total  area 

C4SWEEP  is  the  quarter  chord  sweep  of  the  wing 

ZCFC4W  is  the  distance  parallel  to  the  Z  axis,  from  wing 

root  quarter  chord  point  to  fuselage  centerline. 
BMW  is  the  maximum  body  width 
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2.9  TRIM  AERODYNAMICS 

Aircraft  can  be  trimmed  using  any  one  or  a  combination  of 
several  devices.  These  trimming  devices  fall  into  one  of  three 
classes,  engines,  moveable  surfaces  such  as  canards  and  horizontal 
tails,  or  flaps.  The  aerodynamic  methods  for  each  of  these  are  as 
follows . 

2.9.1  Body  Engines  or  Nacelle  Engines 

If  engines  are  equipped  with  thrust  vectoring  devices,  then 
these  devices  can  be  used  to  trim  the  aircraft.  In  this  case  the 
trim  aero  coefficients  are  simply  computed  from  the  thrust  vector 


CLT  =  CT  *  SIN ( ALPHAT ) 

CDT  =  CT  *  COS (ALPHAT) 

CMT  =  CLT  *  XARM/MACW  +  CDT  *  ZARM/MACW 


(2.9.1) 
(2.9.2  ) 
(2.9.3  ) 


where 


ALPHAT 

XAMR 

ZAMR 

MACW 


is  the  lift  coefficient  due  to  trim  on  the  engine 

is  the  drag  coefficient  due  to  trim  on  the  engine 

is  the  pitching  moment  coefficient  due  to  the  trim 
on  the  engine 

is  the  thrust  vector  angle  for  the  trim  in  the 
plane  of  symnetry 

is  the  X  moment  arm  for  the  engine 

is  the  Z  moment  arm  for  the  engine 

is  the  mean  aerodynamic  chord  of  the  wing 


2.9.2  All  Moveable  Surfaces 

For  all  moveable  surfaces,  such  as  canards  and  horizontal 
tails,  trim  can  be  effected  simply  by  changing  the  angle  of  attack 
on  the  device.  In  this  case  lift  and  drag  can  be  computed  from 
the  methods  described  in  section  2.7  and  the  pitching  moment  can 
be  found  from 

CMT  =  (CLT  *  XARM/MACS  +  CDT  *  ZARM/MACS )  *  SS  *  MACS/ 

(SREF  *  MACW)  (2.9.4) 

where 

CLT  is  the  lift  coefficient  due  to  trim  on  the  surface 

CDT  is  the  drag  coefficient  due  to  trim  on  the  surface 

XARM  is  the  X  moment  arm  of  the  surface 

ZARM  is  the  Z  moment  arm  of  the  surface 

MACS  is  the  mean  aerodynamic  chord  of  the  surface 

MACW  is  the  mean  aerodynamic  chord  of  the  wing 

SS  is  the  surface  area 

SREF  is  the  reference  area 
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Methods  for  computing  the  trim  coefficients  for  flap  are  take 
from  AEROX  (20).  The  methods  are  as  follows: 

Flap  effectiveness.  The  nonlinear  flap-effectiveness  ratio 
is  the  average  of  the  flap  effectiveness  ratios  at  zero  angle  of 
attack  and  at  the  prescribed  higher  angle  of  attack. 

FER  =  (1.  +  COS ( ALPHA  +  DELF  ) /COS ( ALPHA )  )  *  SQRT ( FLAPC/MACS  * 
FLAPS/ (.3  *  (SPANS  -  BDMAXS  )  )  )  when  MACH  <  2. 


or 


FER  = 

FLAPS  * 

(  1  A  + 

1.4  *  SIN ( ALPHA  +  DELF)2 

)/ ( SS  *  ( 1 

SIN ( ALPHA ) z ) ) 

when 

MACH  _>  2. 

where 

DELF 

is 

the 

flap 

deflection  angle 

FLAPC 

is 

the 

f  lap 

average  chord 

MACS 

is 

the 

mean 

aerodynamic  chord  of  the 

surface 

FLAPS 

is 

the 

flap 

span 

SPANS 

is 

the 

surface  span 

BDMAXS 

is 

the 

average  body  diameter  alone  the  surface 

Flap-induced  lift  coefficient.  The  incremental  lift  coefficient 
accompanying  trailing-edge  flap  deflection  is 

CLFLAP  =  FER  *  CLALPHA  *  DELF  *  FM 

(2.9.6) 


where 

CLALPHA  is  the  surface  lift  curve  slope  (unflapped) 

FM  accounts  for  variations  over  the  Mach  number  range 

FM  =  .25;  when  MACH  £  .4 

FM  =  .  322  -  MACH/6.;  when  .4  <  MACH  <_1.4 

FM  =  . 08 3 5 /SQRT ( MACH2  -  1.);  when  1.4  <  MACH  <  2.0 

FM  =  1 . ;  when  MACH  >2.0 

Flap-induced  drag  coefficient.  The  incremental  drag  coefficient 
accompanying  deflection  of  the  surface  trailing-edge  flaps  is 
assumed  to  be: 


CDFLAP  =  CLFLAP  *  TAN( ALPHA  +  DELF/2 .  ) 


(2.9.7) 


Flap-induced  pitching  moment  coefficient .  The 
pitching  moment  coefficient  produced  by  the  flap  is 


CMFLAP  =  -CLFAP  *  EFLAP/MACS 


CSSSeSBaSHBK 


incremental 

(2.9.8) 


where 

EFLAP  accounts  for  the  eccentricities  of  the  centroids  of  the 
loading  produced  by  the  flaps 

EFLAP  =  LACS  -  LCG  +  MACS/4;  when  MACH  £  .4 

EFLAP  =  LACS  -  LCG  +  MACS  *(.25  +  (MACH  -  .4)  *  .5  - 

FLAPC/MACS )  ) ;  when  .4  <  MACH  <  1.4 

EFLAP  =  LACS  -  LCG  +  MACS  *  (.75  -  FLAPC/MACS  )  ; 

when  MACH  >  1.4 

where 

LACS  is  the  X  location  of  the  surface  aerodynamic  center 

LCG  is  the  X  location  of  the  aircraft  center  of  gravity 

The  value  of  the  pitching  moment  must  now  be  referenced  to 
the  aircraft.  This  is  accomplished  by  the  factor. 

CMFLAP  =  CMFLAP  *  SS  *  MACS/(SREF  *  MACW )  (2.9.9) 


2.10  PITCHING  MOMENT  CURVE  SLOPE  ( CMALPHA ) 


According  to  Roskam  (23),  CMALPHA  can  be  estimated  for  all 
subsonic  Mach  numbers  by  the  equation 

CMALPHA  =  (XCG  -  XAC )  *  CLALPHA  (2.10.1) 

where 

XCG  is  the  X  location  of  the  aircraft  center  of  gravity 

XAC  is  the  X  location  of  the  aircraft  aerodynamic  center 

CLALPHA  is  the  lift  curve  slope  of  the  aircraft  determined  by 

CLALPHA  =  CLALPWB  +  CLALPC  +  CLALPH  *  (1  -  DWF ) 

where  the  lift  curve  slopes  for  wing-body,  canard,  and  nacelles 
are  given  by 


CLALPWB  =  KWB  *  CLAW  *  TNW  *  SW/SREF 
CLALPC  =  CLAC  *  TNC  *  SC/SREF 
CLALPH  =  CLAH  *  TNH  *  SH/SREF 


KWB  =  1 

-  .25  » 

*  ( BDW/ SPANW ) 2  +  .025 

*  BDW/ SPANW 

and 

BDW 

is 

the 

body  width 

SPANW 

is 

the 

wing  span 

CLAW 

is 

the 

wing  lift  curve 

slope  from  section  2.7 

CLAC 

is 

the 

canard  lift  curve  slope  from  section  2.7 

CLAH 

is 

the 

horizontal  tail 

lift  curve  slope 

from  section  2 

sw 

is 

the 

wing  area 

SC 

is 

the 

canard  area 

SH 

is 

the 

horizontal  tail 

area 

TNW 

is 

the 

wing  dynamic  pressure  ratio  from 

section  2.8 

TNC 

is 

the 

canard  dynamic  pressure  ratio  from  section  2.8 

TNH 

is 

the 

horizontal  tail 

dynamic  pressure 

ratio  from 

section  2 . 8 

SREF 

is 

the 

reference  area 

DWF 

is 

the 

downwash  factor 

from  the  wing  to 

the  horizontal 

tail 


2.11  THE  VARIATION  OF  PITCHING  MOMENT  COEFFICIENT  WITH  ANGLE  OF 
ATTACK  RATE  ( CMALPD ) 

Roskam  (23)  states  that  CMALPD  can  be  expressed  as  the  sum  of 
components  for  canards  and  horizontal  tails  because  wing  and  body 
contributions  are  small.  With  this  CMALPD  can  be  expressed  as 


CMALPD  =  CMACD  +  CMAHD 
with 

CMACD  =  -2.  *  CLALPC  *  (XCG  -  XACC)2 


CMAHD  =  -2.  *  CLALPH  *  ( XACH  -  XCG)2  *  DWF 
where 


(2.11.1 


(2.11.2) 


(2.11.3) 


CLALPC  is  the  canard  component  lift  curve  slope  from  section 
2.10 

CLALPH  is  the  horizontal  tail  component  lift  curve  slope  from 
section  2.10 

KCG  is  the  X  location  of  aircraft  center  of  gravity 

XACC  is  the  X  location  of  the  canard  aerodynamic  center 

XACH  is  the  X  location  of  the  horizontal  tail  aerodynamic 

center 

DWF  is  the  down  wash  factor  from  the  wing  the  the  horizontal 

tail 
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2.12  THE  VARIATION  OF  PITCHING  MOMENT  COEFFICIENT  WITH  PITCH  RATE 

( CMQ ) 

The  methods  for  calculation  of  CMQ  are  taken  from  Roskam 
(23).  The  derivative  CMQ  may  be  considered  to  be  the  sum  of  a 
wing,  canard,  and  a  tail  contribution,  the  effect  of  the  fuselage 
being  usually  small: 

CMQ  =  CMQW  +  CMQC  +  CMQH  (2.12.1) 

For  the  canard  contribution: 

CMQC  =  -2.  *  CLALPC  *  ( XCG  -  XACC)2  (2.12.2) 

where 

CLALPH  is  the  canard  component  of  lift  curve  slope  from 
section  2.10 

XCG  is  the  X  location  of  the  aircraft  center  of  gravity 

XACC  is  the  X  location  of  the  canard  aerodynamic  center 

For  the  horizontal  tail  contribution: 

CMQH  =  -2  *  CLALPH  *  ( XACH  -  XCG)2  (2.12.3) 

where 

CLALPC  is  the  horizontal  tail  component  of  the  lift  curve 
slope  from  section  2.10 

XACH  is  the  X  location  of  the  horizontal  tail  aerodynamic 
center 

For  the  wing  contribution 

CMQW  =  CMQWMO  *  (R1/R2  +  3./B)/(R3/R4  +3.)  (2.12.4) 

with 

CMQWMO  =  CLAW  *  C0S(C4SWEEP)  *  (T1/T2  +  .041667  *  T3/T4  +  .125) 
and 

B  =  SQRT( 1 .  -  MACH2  *  C4SWEEP ) 

R 1  =  ARW3  *  TAN ( C 4 SWEEP ) 2 

R2  =  ARW  *  B  +  6.  *  COS ( C4SWEEP ) 

R3  =  R 1 

R4  =  ARW  +  6.  *  COS ( C4SWEEP ) 

T 1  =  ARW  *  (2.  *  XW2  +  .5  *  XW) 

T2  =  ARW  +  2.  *  C0S( C4SWEEP ) 

T3  =  ARW3  *  TAN ( C4SWEEP  ) 2 
T4  =  ARW  +  6.  *  COS ( C4SWEEP ) 

XW  =  XACW  -  XCG 


CLAW 

is 

C4SWEEP 

is 

ARW 

is 

XACW 

is 

XCG 

is 

the  lift  curve  slope  of  the  wing 
the  quarter  chord  sweep  of  the  wing 
the  aspect  ratio  of  the  wing 

the  X  location  of  the  aerodynamic  center  of  the  wing 
the  X  location  of  the  aircraft  center  of  gravity  ' 
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2.13  HIGH-LIFT  SYSTEM  AERODYNAMICS 

The  empirical  method  for  predicting  lift,  drag,  and  moment  of 
an  airplane  with  flaps  and  slats  deployed  consist  of  adding  the 
incremental  effects  of  the  high  lift  system  to  the  clean  airplane 
aerodynamics.  Figure  2.13.1  illustrates  the  manner  in  which  the 
incremental  effects  of  a  flap  can  be  applied  to  the  clean-wing 
aerodynamics.  The  following  techniques  for  estimating  these 
increments  were  taken  from  the  Large  Aircraft  Program  ( 1 ) . 

2.13.1  Aerodynamic  Characteristics  of  Two-Dimensional  Hicrh-Lift 


Devices 

The  lift  effectiveness  of  plain  trailing-edge  flaps  can  be 
estimated  from  thin-airfoil  theory.  The  rate  of  change  of  lift 
with  flap  deflection  at  a  constant  angle  of  attack  is  given  by 


cio  =  2  *  (of  +  sin(of; 


(2.13.1 ) 


with 

COS (OF)  =  1  -  2  *  CFP 
where 

CPF  is  the  flap  to  wing  chord  ratio 

This  equation  is  plotted  in  figure  2.13.1  as  a  function  of 
flap-chord  ratio.  The  theory  considers  only  a  bent  flat  plate  and 
does  not  include  effects  of  thickness  or  large  deflection  angles. 
The  effects  are  accounted  for  in  Reference  34  and  35  by  empirical 
flap  efficiency  factors,  as  shown  in  figures  2.13.2  through 
2.13.5.  The  lift  of  a  plain  flap  may  now  be  expressed  as 


DCL  =  ETA  *  CIO  *  DFR 


where 


(2.13.2) 


ETA  is  the  plain-flap  efficiency  factor  from  figure  2.13.2 
depending  on  the  flap  deflection  angle  DFR,  plus  the 
included  angle  of  the  flap  trailing  edge  PHI. 

CIO  is  the  rate  of  change  of  lift  with  flap  deflection  at 
constant  angle  of  attack  from  equation  2.13.1 
DFR  is  the  flap  deflection  angle  in  radians 

This  procedure  is  extended  to  slotted  flaps  with  Fowler 
motion  by  evaluating  CIO  at  a  flap-chord  ratio  based  on  the 
extended  chord.  For  double  -  or  triple  -  slotted  Fowler  flat 
segments  the  lift  increment  is  obtained  by  summing  the  incremental 
lift  increments  for  each  flap  segment.  The  result  is 


/"f  ETA(i)  *  C 1  0 (  i  )  *  DFR(i) 
i=  1 


2.13.3) 


DCF1 
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where 

i  is  a  subscript  that  indicates  the  1st,  2nd,  and  3rd 

flap  segment  of  the  slotted  flap 
I  is  the  number  of  slots  or  segments  in  the  flap  segment 

ETA(i)  is  the  slotted-flap  efficiency  factor  from  Figures 
2.13.3,  2.13.4,  or  2.13.5  for  the  ith  flap  segment 
C10(i)  is  the  lift  for  the  ith  flap  segment 
DFR  is  the  flap  deflection  of  the  ith  flap  segment 


Fi jure  2.13.-5  Turning  Efficiency  of  Plain  Trailing-Edge  Flaps 
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The  method  of  summation  and  the  geometry  definition  required 
to  evaluate  equation  2.13.3  is  shown  in  figure  2.13.6. 

The  effects  of  leading-edge  high-lift  devices  on  the  wing  lift  at 
zero  angle  of  attack  is  estimated  from  thin-airfoil  theory  as 

C1LE  =  2  *  ( SIN ( OS )  -  OS)  (2.13.4) 

where 

COS (OS)  =  1  -  2  *  CSOS 
and 

CSOS  is  the  slat  to  wing  chord  ratio 

Unlike  trailing-edge  flaps,  the  deflection  of  a  nose  flap  causes  a 
loss  in  lift  at  zero  angle  of  attack.  The  increment  in  lift,  DC1S 
is 

DC  IS  =  C1LE  *  DSR  (2.13.5) 

where  DSR  is  the  leading  edge  flap  deflection  angle  in  radians, 
positive  nose  down. 

The  two-dimensional  maximum  lift  increment,  DCLM,  due  to  a 
trailing  edge  plain  flap  deflection  is  given  in  reference  34  as 

DCLM  =  Kt  *  Kd  *  DCL  *  RMAX  (2.13.6) 

and  similarly,  for  single,  double  -  or  triple  -  slotted  flaps  the 
equation  is 

I 

DC1MF  =  Kt  *  Kd  *  £?  DCL(i)  *  RMAX(i)  (2.13.7) 

i=  1 

where 

DCL(i)  is  the  predicted  lift  increment  determined  for  the  ith 
flap  segment  from  equation  2.13.7 
Kt  is  taken  from  the  plot  in  figure  2.13.7 

Kd  is  taken  from  the  plot  in  figure  2.13.8 

RMAX  is  the  theoretical  relationships  between  DCLM  and  DCLO 

given  in  reference  35  as 

=  1.  -  (OF/ (OF  +  SIN(OF)))  *  (1.  +  ALOG ( . 5  *  (X  +  OF))/ 

SIN ( . 5  *  X-OF ))))/( OF  *  TAN ( . 5  *  X)))  (2.13.8) 


RMAX 


LEADING  EDGE  RADIUS /THICKNESS  RATIO 


Figure  2.1 


Fi aure  2.13 


3.7  Maximum  Lift  Correlation  Factor  for  Trailing 
Edge  Flaps 


. 8  Flap  Angle  Correlation  Factor  Versus  Flap 
De  f lec  t ion 


COS(X)  =  2.0  *  CSOS 


-  1.0 

and 

COS ( OF )  =  1.0  -  2.0  *  CFP 
with 

XSOC  is  the  position  of  separation 

CFP  is  the  sum  of  the  flap  to  wing  chord  ratios 

The  equation  relates  the  theoretical  maximum-lift  increment 
to  the  chord  of  the  flap  and  the  position  of  separation,  XSOC,  on 
the  airfoil.  The  choice  of  the  separation  point,  XSOC,  to 
determine  the  maximum-lift  ratio  from  equation  2.13.8  depends  on 
the  leading-edge  configuration.  For  clean  leading-edge  airfoils, 
the  point  of  flow  separation  is  assumed  at  the  leading  edge, 

XSOC  =  0.  For  airfoils  with  leading-edge  high-lift  devices,  the 
point  of  flow  separation  is  assumed  to  be  at  the  knee  of  the 
leading  edge  device,  XSOC  =  CSOS. 

The  two-dimensional  increment  in  maximum-lift  coefficient  of 
leading-edge  devices  is  predicted  in  reference  34  as  follows: 

DC  IMS  =  CMXLE  *  EMAX  *  ED  *  DSR  (2.13.9) 

where  according  to  thin-airfoil  theory, 

CMXLE  =  2.0  *  SIN(OS) 

COS(OS)  =  ( 1 .  -  2.  *  CSOC) 

CMXLE  is  presented  in  figure  2.13.9 
EMAX  is  taken  from  the  plot  in  figure  2.13.10 

ED  is  taken  from  the  plot  in  figure  2.13.11 

DSR  is  the  slat  deflection 

The  maximum-lift  efficiency  factor,  EMAX,  depends  on  the  type 
of  leading-edge  device  and  on  the  ratio  of  leading-edge  radius  to 
maximum  airfoil  thickness;  ED  is  an  efficiency  factor  that 
accounts  for  large  leading-edge-flap  deflections. 

In  the  case  of  two-dimensional  moment  increment,  the 
methodology  for  predicting  the  pitching  moment  is  developed 
parallel  to  the  methods  used  for  estimating  the  lift  increment, 
which  extends  thin-airfoil  theory  to  cover  multiple-slotted  flaps 
with  extendable  chords.  The  trai ling-edge-f lap  pitching-moment 
increment  at  zero  angle  of  attack  is  given  in  reference  34  as: 

DCM  =  DCL  *  RCM  *  Km  (2.13.10) 


85 


."'■.'•VTVFT i  ~~*.w  >L’  *.* 


DCL  is  the  predicted  lift  increment  for  either  trailing- 

edge  or  leading-edge  devices 

RCM  is  the  theoretical  center-of -pressure  location  from  thin- 

airfoil  theory  (figure  2.13.12) 

Km  is  an  empirical  factor  developed  from  experimental  data 

( figure  2.13.13) 

In  the  case  of  the  profile  drag  increment,  flap  drag 
increments  at  ALPHA  =  0  for  plain  and  single-slotted  flaps  are 
obtained  from  Figure  2.13.14  and  2.13.15.  These  figures  were 
obtained  from  Section  6.1.7  in  the  DATCOM.  For  double-  and 
triple-slotted  flaps,  figure  2.13.16  is  used  to  obtain  the 
alpha  =  o  drag  increment. 
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Figure  2.13.  i.:  Flap  Center  of  Pressure  Location  as  given 
by  Thin  Airfoil  Theory 
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Figure  2.13.11*  rwo-Dimensional  Drag  Increment  Due  to  Plain 
Flaps 


2.13.2  Lift  of  High-Lift  Devices 


The  untrimmed  equation  for  lift  can  be  expressed  as 

CL  =  CLO  +  CLALPHA  *  ALPHA  (2.13.11) 

The  increment  in  CLO  and  CLMAX  caused  by  a  trailing-edge- 
flap  deflection  can  be  estimated  from 

DCLOF  =  DC1F  *  CL0C1  *  Kc  *  Kb  (2.13.12) 


DCLMF  =DC 1 MF  *  CL001  *  Kc  *  Kb  *  COS(SWEEPMT)  (2.13.13) 

where 

DCIF  is  defined  in  section  2.13.1 

DCIMF  is  defined  in  section  2.13.1 

Kc  is  taken  from  the  plot  in  figure  2.13.18 

Kb  is  taken  from  the  plot  in  figure  2.13.19 

The  factor  CL0C1 ,  determined  from  the  Polhamus  lift  equation 

CL0C1  =  ARW/ (2.0  +  SQRT( 4 . 0  +  ( SW  *  ARW/SPLANX/COS ( SWEEPMT ) ) 2 ) ) 

(2.13.14) 

where 

ARW  is  the  aspect  ratio  of  the  wing 

SW  is  the  area  of  the  wing 

SPLANX  is  the  planform  area  of  the  wing  with  flaps  extended 

SWEEPMT  is  the  maximum  thickness  sweep  of  the  wing 

and  converts  the  sectional  values  DCIF  and  DCIMF  to  the  three 
dimensional  case.  The  sectional  values  are  either  obtained  from 
input  to  the  program  or  generated  internally  for  certain  types  of 
high-lift  system.  The  method  of  generating  the  section  values  is 
discussed  in  section  2.13.1. 
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The  increment  in  lift  at  zero  angle  of  attack  is 
approximately  zero  when  a  slat  is  deflected.  The  slat  acts  to 
delay  separation  from  the  wing  leading  edge  and  thus  allows  higher 
angles  of  attack  and,  consequently,  higher  values  of  maximum  lift 
before  the  wing  stalls.  An  estimate  of  the  increase  in  maximum 
lift  of  a  slat  is  represented  by 


DCLMS  =  DC  IMS  *  CL0C1  *  Ks  *  COS (SWEEPLE) 


(2.13.15) 


where 

SWEEPLE  is  the  leading  edge  sweep  of  the  wing 

Ks  is  taken  from  the  plot  in  figure  2.13.20 

DC  1  MS  is  the  sectional  increase  in  maximum  lift  of  a  slat 

2.13.3  Drag  of  High-Lift  Devices 


The  untrimmed  equation  for  drag  can  be  expressed  by 
CD  =  CDO  +  K  *  ( CL-DCL ) 2 


(2.13.16) 


where  the  drag-due-to  lift  factor,  K,  with  high  lift  devices  can 
be  estimated  by 


K  =  ( CD-CD0/CL2 )  *  ( SREF/SPLANX ) 


(2.13.17 


The  change  in  minimum  drag  for  the  high-lift  configuration  is 
expressed  as 


DCDMIN  =  DCDF  +  DCDS  +  CDI  +  CDCLEAN 


(2.13.18) 


where  increments  are  summed  for  the  profile  drag  caused  by  the 
flaps  (DCDF)  and  the  slats  (DCDS).  The  increment  in  profile  drag 
of  the  flaps  and  slats  can  be  estimated  from  sectional  drag  data, 
using 


DCDF  =  DCDFS  *  COS(SWEEPHL)  *  Kd 


(2.13.19) 


DCDS  =  DCDSS  *  COS (SWEEPLE)  *  Kd 


(2.13.20) 


where 


SWEEPHL  is  the  sweep  of  the  high-lift  system 

Kd  is  taken  from  the  plot  in  figure  2.13.21 

DCDFS  is  the  section  DCDF  given  in  section  2.13.1 

DCDSS  is  the  section  DCDS  given  in  section  2.13.1 

Deflection  of  a  flap  procedures  an  increase  in  lift  at  zero 
angle  of  attack  which  in  turn  produces  an  induced  drag  given  by 

CDI  =  Ka  *  Kf  *  DCLOF2/(PI  *  ARW  5 


(2.13.21  ) 


where 

Ka  is  taken  from  the  plot  in  figure  2.13.22 

Kf  is  taken  from  the  plot  in  figure  2.13.22 

Ka  and  Kf  are  factors  which  account  for  the  noneleptical  span 
loading  of  partial-span  flaps. 

The  deflection  of  a  flap  increases  the  camber  of  the  airfoil. 
In  reference  33,  thin-airf oi 1  camber  theory  is  used  to  relate  the 
displacement  of  a  polar  with  flaps  to  the  lift  increment  of  the 
flap  at  zero  angle  of  attack  by  the  equation. 

DCL  =  DCLCLEAN  +  DCLF 

where 

DCLF  =  (1  -  1 / ( PI  *  ARW  *  AK))  *  DCLOF/ ( 1  +  1.16  *  CLOCl  * 

( .5  -  CF)  (2.13.22) 

and 

AK  is  the  drag  due  to  lift 

CF  is  the  flap  to  wing  chord  ratio 

DCLOF  is  the  increment  in  CLO  caused  by  a  flap  deflection  from 

Section  2.13.2 


2.13.4  Moment  of  High-Lift  Devices 


The  pitching  moment  increment  caused  by  a  flap  on  a  swept 
wing  is  represented  by 

DCMAF  =  DCMFS  *  CPF2  *  Km  +  . 5  *  ARW  *  TAN(SWEEPC2)  *  DCIF  *  Ksw 

(2.13.23) 


where 


DCMFS  is  the  sectional  change  in  moment  at  zero  alpha  due  to 

flap  deflection 

DCIF  is  the  sectional  change  in  lift  at  zero  alpha  due  to 

flap  deflection 

CPF  is  the  flap  to  wing  chord  ratio 

ARW  is  the  aspect  ratio  of  the  wing 

SWEEPC2  is  the  half  chord  sweep  of  the  wing 
Km  is  taken  from  the  plot  in  figure  2.13.23 

Ksw  is  taken  from  the  plot  in  figure  2.13.24 
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3.0  STABILITY  COEFFICIENTS  METHODS 

The  methods  described  below  are  used  to  compute  the  lateral 
stability  derivatives  of  an  aircraft.  The  methods  include  some 
taken  from  DATCOM  (8),  and  some  that  have  been  developed  for  use 
in  SACP.  The  derivatives  that  are  to  be  computed  are;  the 
variation  of  rolling  moment  coefficient  with  roll  rate  (CLP),  the 
variation  of  rolling  moment  coefficient  with  aileron  angle  (CLDA), 
the  variation  of  rolling  moment  coefficient  with  sideslip  angle 
( CLB ) ,  the  variation  of  yawing  moment  coefficient  with  sideslip 
angle  ( CNB ) ,  the  variation  of  yawing  moment  coefficient  with 
rudder  angle  (CNDR),  and  the  variation  of  yawing  moment 
coefficient  with  yaw  rate  (CNR). 

3.1  CLP 

The  variation  of  rolling  moment  coefficient  with  roll  rate 
(CLP)  can  be  computed  as  the  simple  sum  of;  CLP  for  the  Body  with 
wing  interference  (CLPWB),  CLP  for  the  horizontal  tail  ( CLPHT ) , 

CLP  for  the  canard  ( CLPC )  and  CLP  for  the  vertical  tail  ( CLPVT ) . 
The  equation  for  CLP  is  then  just 

CLP  =  CLPWB  +  CLPW  +  CLPHT  +  CLPC  +  CLPVT  (3.1.1) 

where  computation  of  each  of  the  components  is  described  below. 

3.1.1  CLP  For  Lifting  Surfaces 

The  method  for  computing  the  wing  rolling  derivative,  CLPW, 
is  given  in  the  DATCOM.  This  method  can  be  extended  to  other 
lifting  surfaces  by  referencing  the  results  to  the  appropriate 
reference  area.  This  section  presents  the  DATCOM  methods  for 
estimating  the  CLPS  at  subsonic  and  supersonic  speeds  for  all 
lifting  surfaces. 

The  DATCOM  method  for  computation  of  CLPS  at  subsonic  speeds 
accounts  for  the  variations  in  wing  lift  curve  slope,  drag  due  to 
lift,  and  profile  drag,  as  well  as  the  effect  of  dihedral.  The 
values  of  CLPS  at  a  given  lift  coefficient  at  subsonic  speeds, 
based  on  the  product  of  the  surface  area  and  the  square  of  the 
surface  span,  SA  *  SPAN2,  is  given  by 

CLPS  =  ( RDPZL  *  (K/BETA)  *  ( CLALP/CLALPO )  *  DIHP)  +  DCLPD  (3.1.2) 
where 

RDPZL  is  the  damping  parameter  at  zero  lift  obtained  from 

figure  3.1.1  as  a  function  of  SWPBETA  and  BA/K 
SWPBETA  is  the  compressible  sweep  parameter  given  by 


SWPBETA  =  ARCTAN( TAN( C4SWEEP ) /BETA ) 


SUBSONIC  SPEEDS 


(b)  X  =  0.25 

(3A 


K 


Fioure  3.1.1 
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BETA  =  SQRTM.O  -  MACH2) 


BAK  =  BETA  *  ARS/K 

ARS  is  the  aspect  ratio  of  the  surface 

The  parameter  K  is  the  ratio  of  the  two-dimensional  lift  curve 
slope  at  the  appropriate  Mach  number  to  2.0  *  PI/BETA.  i.e., 

K  =  CLALP/ (2.0  *  PI/BETA)  (3.1.3) 

CLALP  is  the  surface  lift  curve  slope  at  any  lift  coefficient 
below  stall 

CLALPO  is  the  surface  lift  curve  slope  at  zero  lift 
DIHP  is  the  dihedral-effect  parameter  given  by 

DIHP  =  1.0  -  2.0  *  (ZCG/( SPAN/2.0) )  *  SIN(DIHSR) 

+  ((  ZCG/( SPAN/2.0)2  *  ( ( SIN(DIHSR  )2  )  )  )  (3.1.4) 

DIHSR  is  the  geometric  dihedral  angle,  positive  for  surface 
tip  above  surface  root. 

ZCG  is  the  vertical  distance  between  the  aircraft  center  of 

gravity  and  the  surface  root  chord,  positive  for  c.g. 
above  the  root  chord 

DCLPD  is  the  increment  in  the  roll  damping  derivative  due  to 
drag.  It  is  given  by 

DCLPD  =  DDLRDP  *  CL2  -  1 /8  *  CDO  (3.1.5) 

where 

DDLRDP  is  the  drag  due  to  lift  roll  damping  parameter  obtained 
from  figure  3.1.2  as  a  function  of  ARS  and  C4SWEEP 
CL  is  the  surface  lift  coefficient 

and 

CDO  is  the  total  zero  lift  drag  coefficient 

For  wings  of  low  aspect  ratio  and/or  high  sweep  the  accuracy 
of  the  method  rapidly  deteriorates  with  increasing  CL,  even  when 
experimental  values  of  lift  and  drag  are  used.  The  error  results 
from  the  fact  that  the  high  values  obtained  from  figure  3.1.2  for 
these  configurations  are  not  realized  in  practice.  Therefore,  as 
CL  increases  the  calculated  values  of  the  roll  damping  derivative 
become  progressively  smaller  than  those  given  by  experiment. 

At  supersonic  speeds,  design  charts  based  on  theoretical 
calculations  are  presented  for  estimating  CLP  of  surfaces.  In 
this  case  CLP  is  given  by 


1  OR 


CLP  =  CLPATF  *  ARS  *  CLPCLT 


(3.1.6) 


SUBSONIC  SPEEDS 


where 

CLPAT  is  the  theoretical  roll-damping  parameter  obtained  from 

figure  3.1.3a  through  3.1.3e 

and 

CLPCLPT  is  the  empirical  thickness  correction  factor  obtained 
from  f igure  3.1.4 

The  sonic  trailing  edge  boundaries  on  figure  3.1.3a  through 
3.1.3e  represent  an  upper  limit  for  the  true  theoretical  values  of 
the  derivatives.  Values  below  the  sonic  trailing-edge  boundary 
are  for  wings  with  subsonic  trailing  edges  and  are  in  violation  of 
one  of  the  basic  assumptions  of  the  theory.  For  configurations 
with  subsonic  trailing  edges  the  design  charts  wil]  over-estimate 
the  roll  damping  derivative. 

It  should  be  noted  that  the  "kinks"  in  the  curves  of  figure 
3.1.3a  through  3.1.3e  correspond  to  the  conditions  of  sonic 
leading  edges.  Experimental  evidence  shows  that  these  kinks  do 
not  occur  in  practice. 

3.1.2  CLP  For  the  Vertical  Tail 


The  contribution  of  CLP  due  to  vertical  tails  is,  according 
to  Roskam  ( 23  )  , 

CLPVT  =  -2.0  *  CLALPVO  *  (ZV/SPAN)2  *  NV  *  SV/SW  *  VNO  (3.1.7) 
where 


CLALPVO  is  the  vertical  tail  lift  curve  slope 

ZV  is  the  vertical  location  of  the  mean  aerodynamic  chord 

of  the  vertical  tail 

NV  is  the  vertical  tail  sidewash  factor 

SV  is  the  area  of  the  vertical  tail 

SW  is  the  area  of  the  wing 

SPAN  is  the  vertical  tail  span 

VNO  is  the  number  of  vertical  tails 


3.1.3  CLP  for  Wing-Body  Interaction 


This  section  presents  methods  for  estimating  the  wing-body 
contribution  to  the  rolling  derivative  CLP  at  subsonic  and 
supersonic  speeds.  This  derivative  is  the  change  in  rolling- 
moment  coefficient  with  change  in  the  wing  tip  helix  angle  and  is 
referred  to  as  the  roll  damping  derivative. 


For  subsonic  speeds  there  is  no  body  contribution  to  CLPWB  so 
the  value  is  simply 


1  1  0 


-  1 1 


SUPERSONIC  SPEEDS 


(CONTDl 


wws 


CLPWB  =  CLPW 


CLPW  is  the  same  as  calculated  in  section  3.1.1 

The  CLPWB  for  supersonic  conditions  is  given  in  the  DATCOM  as 


CLPWB  =  CLPW  *  SW/SREF  *  BCF 


where 


(3.1.8) 


BCF  is  the  semi-empirical  body-correction  factor  obtained 

from  figure  3.1.5  as  a  function  of  the  ratio  of  maximum 
body  diameter  to  wing  span  and  the  Mach  number  normal  to 
the  leading  edge. 
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3.2  CLDA 

The  variation  of  rolling  moment  coefficient  with  aileron 
deflection,  CLDA,  can  be  computed  by  a  procedure  in  Roskam  (23)  as 

follows . 

STEP  1 

Obtain  the  rolling-moment  effectiveness  parameter  BCLPDK  of 
two  full-chord  controls  (Cf/C  =  1.0)  anti-symmetr ically  deflected, 
as  a  function  of  BA/K  and  SWPBETA  from  Figure  3.2.1. 

The  parameter  K  is  the  ratio  of  the  two-dimensional  lift- 
curve  slope  at  the  appropriate  Mach  number  (CLALP)  to  2  *  PI/BETA, 
i.e.  CLALPM/ ( 2  *  PI/BETA).  For  wings  with  airfoil  sections 
varying  in  a  reasonably  linear  manner  at  the  Mach  of  the  flapped 
portion  of  the  wing  is  adequate.  The  parameter  SWPBETA  is  the 
compressible  sweep  parameter,  given  as: 

SWPBETA  =  ATAN ( TAN ( C4 SWEEP ) /BETA )  (3.2.1) 

Figure  3.2.1  gives  directly  the  effectiveness  parameter  for 
control  spans  measured  from  the  plane  of  symmetry  outboard.  For 
partial  span  controls  having  the  inboard  edge  of  the  control  at 
spanwise  station  N(i)  and  the  outboard  edge  at  N(o)  the 
effectiveness  parameter  is  obtained  as  illustrated  in  Figure 


STEP  2 

Determine  the  rolling  effectiveness  of  two- f u 1 1-chord 
controls  anti-symmetr ically  deflected  by: 

CPLD  =  K/BETA  *  (BCLPDK)  (3.2.2) 

STEP  3 

Determine  the  rolling  effectiveness  of  the  par tia 1 -chord 
controls  (Cf/C  =  constant)  anti-symmetr ica 1 ly  deflected  by: 

CLD  =  ABAD  *  CPLD  (3.2.3) 

where 

ABAD  is  the  absolute  value  of  the  section  lift  effectiveness, 

AD.  AD  is  obtained  from  equation  3.2.4  for  the  particular 
control  under  consideration.  For  anti-symmetric  control 
deflections,  the  value  of  AD  is  based  on  the  deflection 
of  the  surface 


ABAD  =  CLD/CLALPH 


(3.2.4) 


I 


where 


CLD  =  CLDCLDT  *  CLDT  (3.2.5) 

found  in  figures  3.2.4  and  3.2.3  respectively  and  where 

CLALPH  is  the  average  section  lift  curve  slope  over  the  aileron 
span  of  the  wing 

STEP  4 

The  effect  of  a  differential  control  deflection  is  taken  into 
account  by  considering  CLD  of  each  control  as  one-half  the  anti- 
symmetrical  value  (equation  3.2.3)  where  ABAD  is  considered 
separately  for  each  control  and  based  on  its  respective 
def lection . 

Then  the  total  rolling-moment  coefficient  for  differential- 
control  deflection  is  obtained  by: 

CL  =  ( CLD ( L  )  / 2 . 0  *  AIL(L)  -  CLD(R)/2.0  *  AIL ( R  )  (3.2.6) 

where  it  should  be  kept  in  mind  that  a  positive  control  deflection 
is  trailing  edge  down,  and  AIL(L)  and  AIL(R)  are  the  left  and 
right  aileron  deflections.  "The"  aileron  deflection  is  defined  as 

DA  =  1/2(AIL(L)  -  AIL ( R  )  )  (3.2.7) 

it  follows  that 


CLDA  =  CLD(L)  +  CLD(R) 


(3.2.8) 
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3 . 3  CLB 

The  variation  of  rolling  moment  coefficient  with  sideslip 
angle,  CLB,  is  shown  in  Roskam  (23)  to  be  expressed  as  the  sum  of 
the  component  CLB’s  as: 


CLB  =  CLBWB  +  CLBC  +  CLBH  +  CLBV 


(3.3.1  ) 


where 

CLBWB  is  the  CLB  for  wing-body  interaction 

CLBC  can  be  treated  as  the  canard-body  interaction  CLB 

CLBH  can  be  treated  as  the  horizontal  tail-body  interaction  CLB 

CLBV  is  the  vertical  tail  CLB 


CLB  for  all  of  the  lifting  surfaces  can  be  computed  by  the  same 
method,  while  CLBV  must  be  treated  separately. 

3.3.1  CLB  For  Vertical  Tails 

According  to  Roskam  (23)  CLB  for  the  vertical  tail  can  be 
expressed  as 


(3.3.2) 


CLBV  =  -CLALPV  *  SWFV  *  (SV/SW)  *  BV  *  VNO  ( 

where 

BV  =  (ZV  *  COS ( ALPHA )  -  LV  *  COS ( ALPHA )) /SPANW 

ZV  is  the  vertical  location  of  the  vertical  tail  mean 

aerodynamic  chord 

LV  is  the  axial  location  of  the  vertical  tail  mean 

aerodynamic  chord 

CLALPV  is  the  vertical  tail  zero  alpha  lift  curve  slope 

SWFV  is  the  sidewash  factor  of  the  vertical  tail 

SV  is  the  area  of  the  vertical  tail 

SW  is  the  area  of  the  wing 


VNO  the  number  of  vertical  tails 

3.3.2  CLB  For  Lifting  Surfaces 

The  subsonic  rolling  moment  due  to  sideslip,  based  on  the 
product  of  the  surface  span  and  area,  for  straight-tapered 
surface-body  combinations  at  low  angles  of  attack  is  given  by  the 
DATCOM  as; 


(  ( CLSB  *  SI  )  +  ( DIHSD  *  S2  )  -t  DELTCLB  +  S3)  *  57.3 


(3.3.3) 
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51  =  CLBCLC2  *  KMS  *  FLEF  +  ARCF 

52  =  SDEF  *  KMD  +  CLBDIH 

53  =  TWISTS  *  TAN ( C4SWEEP )  *  DELCLBH 

where 

CLSB  is  the  lift  coefficient  of  the  body  and  surface 

CLBCLC2  is  the  surface  sweep  contribution  obtained  from  figure 

3.3-1 

KMS  is  the  compressibility  correction  to  the  sweep 

contribution  obtained  from  figure  3.3.2 
FLEF  is  the  empirical  fuselage-length-effect  correction 

factor  obtained  from  figure  3.3.3 
ARCF  is  the  aspect  ratio  contribution  obtained  from  figure 

3.3.4 

DIHSD  is  the  geometric  dihedral  angle  in  degrees,  positive 

for  the  surface  tip  above  the  plane  of  the  root  chord 
SDEF  is  the  surface  dihedral  effect  for  uniform  geometric 

dihedral  obtained  from  figure  3.3.5 
KMD  is  the  compressibility  correction  factor  to  the  dihedral 

effect  obtained  from  figure  3.3.6 

CLBDIH  is  the  body-induced  effect  on  wing  height  due  to  uniform 
geometric  dihedral  given  by 

CLBDIH  =  -.0005  *  SQRT(ARS)  *  ( DIAEQ/SPANS ) 2 

ARS  is  the  aspect  ratio  of  the  surface 

DIAEQ  is  the  average  body  diameter  at  the  surface  for  non¬ 

body  of  revolution  configuration,  the  average  equivalent 
diameter  at  the  surface  root  should  be  used 

DIAEQ  =  SQRTtaverage  cross  sectional  area/.  07854) 

SPANS  is  the  span  of  the  surface 

DELTCLB  is  the  increment  in  CLB  due  to  the  body-induced  effect 
on  wing  height  for  configurations  with  surfaces  located 
above  or  below  the  midfuselage  height.  The  increment 
given  by 

DELTCLB  =  (1.2  *  ( SQRT ( ARS  ) / 5 7 . 3  )  *  ZCG/SPANS)  *  2  *  DIAEQ/SPANS 
where 

ZCG  is  the  vertical  distance  from  the  center  of  gravity 

to  the  quarter  chord  point  of  the  surface 
TWISTS  is  the  surface  twist  between  the  root  and  tip  sections 
in  degrees  (see  Figure  3.3.7) 

C4SWEEP  is  the  sweep  angle  of  the  quarter  chord  point  of  the 
surface 
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Figure  3.3.5  EFFECT  OF  UNIFORM  GEOMETRIC  DIHEDRAL  ON  WING 
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DELCLBH  is  the  wing  twist  correction  factor  obtained  from 
figure  3.3.7 


The  supersonic  rolling  moment  due  to  sideslip  for  surface- 
body  configurations  with  straight-tapered  surfaces  at  low  angles 
of  attack  is  approximated  in  the  DATCOM  by 

CLB  =  -.061  *  CN  *  CLALPS  *  Tl  *  T2  *  (T3  +T4 )  +  T5  +  DELTCLB 

(3.3.4) 

with 

CN 
Tl 
T2 
T3 
T4 
T5 

where 

TRS  is  the  taper  ratio  of  the  surface 

LESWEEP  is  the  leading  edge  sweep  angle  of  the  surface 

CLALPS  is  the  zero  alpha  lift  curve  slope  of  the  surface 

SDEFS  is  the  surface  dihedral  effect  given  by 

CLPS  is  the  CLP  for  the  surface  from  Section  3.1.1 

and 

SDEFS  =  2/57. 32  *  (1.  0  +  (  2.0  *  TRS )  ) / ( 1  . 0  +  (3.0  *  TRS))  *  CLPS 

The  value  of  CLB  must  now  be  adjusted  for  each  surface  by 

CLBS  =  CLB  *  NS  *  SS/SREF  *  SPANS/SPANW  (3.3.5) 

where 

NS 
SS 

SREF 
SPANS 

and 

SPANW  is  the  span  of  the  wing 


is  the  dynamic  pressure  ratio  of  the  surface 
is  the  reference  area  of  the  surface 
is  the  reference  area  of  the  aircraft 
is  the  span  of  the  surface 


=  CLSB/C0S( ALPHA) 

=  1.0  +  (TRS  *  (1.0  +  LESWEEP)) 

=  (1.0  +  ( LESWEEP/ 2  )  )  *  TAN( LESWEEP) /BETA 
=  MACH2  *  ( COS ( LESWEEP  )  )2/ARS 
=  ( TAN ( LESWEEP ) / 4 ) ^  ‘ 3  2 
=  DIHSD  *  (SDEFS  +  CLBDIH  ) 


3.4  CNB 
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The  variation  of  yawing  moment  coefficient  with  sideslip 
angle,  CNB,  can  be  obtained  from  the  sum  of  the  component  CNB ' s  by 
the  equation 


CNB  =  CNBB  +  CNBN  +  CNBW  +  CNBC  +  CNBH  +  CNBV  +  CNBF  +  CNBHV 

(3.4.1  ) 


where 


CNBB 

CNBN 

CNBW 

CNBC 

CNBH 

CNBV 

CNBF 

CNBHV 


is  the  CNB  of  the  body 
is  the  CNB  of  the  nacelles 


is  the  CNB  of  the  wings 
is  the  CNB  of  the  canards 
is  the  CNB  of  the  horizontal  tails 
is  the  CNB  of  the  vertical  tails 
is  the  CNB  of  the  fins 

is  the  CNB  of  the  horizontal  tail  (vertical  component) 


The  derivatives  for  each  of  these  components  is  described  below. 


3.4.1  CNB  for  Bodies  and  Nacelles 


The  total  body  or  nacelle  sideslip  derivative  CNB  is  given  in 
the  DATCOM  as 


CNB  =  -57.3  *  KN  *  KRL  *  BSA/SREF  *  LB/SPANW 


3.4.2  ) 


where 


BSA 

SREF 

LB 

SPANW 


is  the  empirical  factor  related  to  the  sideslip 
derivative  for  body  plus  wing  interference  obtained 
from  figure  3.4.1  as  a  function  of  geometry 
is  the  empirical  Reynolds-number  factor  obtained  from 
figure  3.4.2 

is  the  side  area  of  the  body  or  nacelle 
is  the  reference  area 
is  the  length  of  the  body  or  nacelle 
is  the  span  of  the  wing 


For  multiple  nacelles,  this  result  is  summed  over  all  of  the 
nacelles . 


3.4.2  CNB  for  Lifting  Surfaces 


The  yawing  moment  of  surface  in  sideslip  is  primarily  caused 
by  the  asymmetrical  induced-drag  distribution  associated  with  the 
asymmetrical  lift  distribution.  The  surface  contribution  to  the 
derivative  CNB  is  important  only  at  large  incidence. 


EMPIRICAL  FACTOR  Kn  RELATED  TO  SIDESLIP  DERIVATIVE  Cn  FOR 
BODY  +  WING-BODY  INTERFERENCE  ** 
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The  DATCOM  method  presented  herein  uses  distinct  equations 
for  computation  of  low  subsonic,  high  subsonic,  and  supersonic 

speeds . 

The  yawing  moment  derivatives  at  low  speeds  is  given  by 

CNBLS  =  CL2  *  ( Q 1  -  TAN ( C 4 SWEEP  )  /Q2 )  *  Q3  (3.4.3) 

with 


Q 1  =  1 . 0/ ( 4 . 0  *  PI  *  ARS ) 

Q2  =  PI  *  ARS  *  (ARS  +  4.0  *  COS ( C4SWEEP )  ) 
Q3  =  COS ( C4SWEEP )  -  ARS/2.0  -  Q4  +  Q5 
Q4  =  ARS2/ (3.0  *  COS ( C4SWEEP  )  ) 

Q5  =  6.0  *  XB  *  S IN ( C4 SWEEP ) / ARS 

where 


CL 

C4SWEEP 

ARS 

XB 


is  the  lift  coefficient  of  the  surface 
is  the  quarter  chord  sweep  of  the  surface 
is  the  aspect  ratio  of  the  surface 

is  the  axial  distance  from  the  surface  aerodynamic 
center  to  the  aircraft  center  of  gravity 


For  subcritical  speeds,  the  low-speed  derivative  can  be 
modified  by  the  Prandt 1 -G lauert  rule  to  yield  approximate 
corrections  for  the  first  order  three-dimensional  effects  of 
compressibility.  The  resulting  expression  is 


CNB  =  CNBLS  *  (Q6/Q7)  *  (Q8/Q9) 


(3.4.4) 


with 


Q6  =  ARS  +  (4.0  *  COS ( C4SWEEP )  ) 

Q7  =  (ARS  *  B )  +  (4.0  *  COS ( C4SWEEP ) ) 

Q8  =  ARS2  *  B2  +  (4.0  *  ARS  *  B  *  COS ( C4 SWEEP )  ) 

-  (8.0  '  COS ( C4SWEEP )  )  2 

Q9  =  ARS2  +  (4.0  *  ARS  *  COS ( C4 SWEEP  )  )  -  (8.0  *  COS ( C4 SWEEP  )  ) 2 
where 

B  =  SQRTd.O  -  MACH2  *  (  COS  (  C4SWEEF  )  )  2 

The  DATCOM  Methods  for  calculation  of  CNB  at  supersonic 
speeds  is  divided  into  methods  for  rectangular  planfcrms,  and 
fully  tapered  sweptback  planforms  with  swept  forward  or  sweptback 
trailing  edges.  The  results  are  mainly  functions  of  planform 
geometry  and  Mach  number.  The  general  trend  of  the  variation  of 
CNB  with  Mach  number  and  aspect  ratio  is  a  reduction  in  the 
magnitude  of  the  derivative  with  an  increase  m  these  parameters. 
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Rectangular  Planforms  ARS  *  BETA  _>  1-0 

The  wing  yawing  moment  due  to  sideslip  for  rectangular  planforms 
referred  to  an  arbitrary  moment  center  is  given  as 

CNB  =  ALPHA2  *  Q 1 0  *  ( Q 1  1  +  Q 1 2  -  PI  *  (Q13  +  Q1 4 )  )  (3.4.5) 

with 

Q10  =  1.0/ (PI  *  ARS2  *  BETA2) 

Q 1 1  =  (4.0  *  MACH2  )  / 3 . 0 
Q 1 2  =  8.0  *  MACH2  *  XOC 1 
Q 1 3  =  (ARS  *  (1.0  ■  BETA2))/BETA 

Q 1 4  =  (3.0  +  BETA2 )/( 3.0  *  BETA2) 

where 

BETA  =  SQRT ( MACH2  -  1.0) 

XOC1  is  the  distance  of  the  origin  of  moments  from  the  mean 
aerodynamic  chord 

Equation  3.4.5  is  valid  for  Mach  number  and  aspect  ratio 
greater  than  that  for  which  the  Mach  line  from  the  leading  edge  of 
the  tip  section  intersects  the  trailing  edge  of  the  opposite  tip 
section  (ARS  *  BETA  _>  1.0) 

Sweptback  Planform 

The  wing  yawing  moment  due  to  sideslip  for  fully  tapered  sweptback 
planforms  is: 

CNB  =  ALPHA2  *  PI/3  *  1 EBECT  *  F9N  +  (ARS2/16  *  FI  IN  +  X0C2  ) 

*  MACH2  *  QBECT )  (3.4.6) 


where 


EBETC 

is 

obtained 

from 

f igure 

3.4.3 

F9N 

is 

obtained 

from 

f igure 

3.4.4 

FI  IN 

is 

obtained 

f  rom 

f igure 

3.4.5 

QBECT 

is 

obtained 

from 

f igure 

3.4.6 

and 

XQC2  is  the  distance  of  the  origin  of  moments  from  the  2/3 
chord  point  of  the  basic  triangular  wing,  measured 
alone  the  longitudinal  axis. 
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3.4.3  CNB  for  Vertical  Tails 

The  variation  of  yawing  moment  coefficient  with  sideslip 
angle  due  to  the  vertical  tail  can  be  represented  as 


CNBV  =  CLALPVO  *  SWFV  *  SV/SREF  *  BV  *  VNO 
where 


(3.4.7) 


CLALPVO  is  the  zero  alpha  lift  curve  slope  of  the  vertical  tail 

SWFV  is  the  side-wash  factor  of  the  vertical  tail 

SV  is  the  area  of  the  vertical  tail 

SREF  is  the  reference  area 

VNO  is  the  number  of  vertical  tails 


BV  = 
LV 

SPANW 

and 

ZV 


(LV  *  COS (ALPHA)  +  ZV  *  S1N( ALPHA )) /SPANW 


(3.4.8) 


is  the  axial  location  of  the  vertical  tail  aerodynamic 
center 

is  the  span  of  the  wing 


is  the  vertical  location  of  the  mean  aerodynamic  chcrd 
of  the  vertical  tail 


3.4.4  CNB  for  Fins 

The  variation  of  yawing  moment  coefficient  with  sideslip 
angle  due  to  fins  is  represented  as 

CNBF  =  CLALPFO  *  SF/SREF  *  BF/SPANW  *  FNO  (3.4.9) 

where 

CLALPFO  is  the  zero  alpha  lift  curve  slope  of  the  fin 

SF  is  the  area  of  the  fin 

SREF  is  the  reference  area 

SPANW  is  the  span  of  the  wing 

FNO  is  the  number  of  fins 

BF  =  LF  *  COS (ALPHA)  +  ZF  *  SIN (ALPHA)  (3.4.10) 

where 


LF 

ZF 


is  the  X  location  of  the  fin  aerodynamic  center 
is  the  Z  location  of  the  fin  mean  aerodynamic  chord 


3.4.5  CNB  for  the  Vert ical  Component  of  the  Horizontal  Tail 

If  the  Horizontal  tails  have  a  slight  dihedral  angle  then 
they  will  have  a  CNB  component  due  to  the  vertical  portion  of  tht 
tail.  This  component  is  defined  as 
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CNBHV 


CLALPH  *  ASEXH/SREF  *  SIN(DIHH)  *  NH  *  BH 


(3.4.11) 


where 

CLALPH  is  the  horizontal  tail  lift  coefficient  due  to  alpha 

ASEXH  is  the  exposed  area  of  the  horizontal  tail 

SREF  is  the  reference  area 

DIHH  is  the  dihedral  angle  of  the  horizontal  tail 

NH  is  the  dynamic  pressure  ratio  at  the  horizontal  tail 

BH  =  ((LH  *  COS ( ALPHA  )  +  ( ZH  *  SIN ( ALPHA )) /SPANW  (3.4.12) 


SPANW  is  the  span  of  the  wing 

LH  is  the  X  location  of  the  aerodynamic  center  of  the 

horizontal  tail 

ZH  is  the  Z  location  of  the  mean  aerodynamic  chord  of  the 

horizontal  tail 
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3 . 5  CNDELTR 


The  variation  of  yawing  moment  coefficient  with  rudder  angle 
CNDELTR  can  be  computed  by 


CNDELTR  =  -CLALPVO  *  SWFV  *  ALFDR  *  SV/SREF  *  BV 
where 


(3.5.1  ) 


ALFDR  is  the  rudder  effectiveness  factor,  obtained  from 

figure  3.5.1 


ALFDR  =  .2923253  +  .4191569  *  CFC  +  .4816274  *  CFG' 

ALFDR  =  1.0  for  an  all  moveable  tail 

ALFDR  =  0.0  for  a  fixed  tail  with  no  rudder 


CFC=  RUDC/MACV 


RUDC 

MACV 

CLALPVO 

SWFV 

SV 

SREF 


is  the  average  rudder  chord 

is  the  mean  aerodynamic  chord  of  the  vertical  tail 
is  the  zero  alpha  lift  curve  slope  of  the  vertical  tail 
is  the  sidewash  factor  of  the  vertical  tail 
is  the  area  of  the  vertical  tail 
is  the  reference  area 


BV  =  (LV  *  COS ( ALPHA )  +  ZV  *  SIN( ALPHA )) /SPANW 


(3.5.2) 


LV 

SPANW 


the  X  location  of  the  vertical  tail  aerodynamic  center 
the  span  of  the  wing 
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3 . 6  CNR 

The  variation  of  yawing  moment  coefficient  with  yaw  rate  CNR 
is  simply  the  sum  of  CNR  for  the  vertical  tail  ( CNRV )  and  CNR  for 
the  wing  body  combination  (CNRWB).  This  is  represented  as 

CNR  =  CNRWB  +  CNRV  (3.6.1) 

where  the  components  are  computed  below. 

3.6.1  CNRWB 


The  variation  of  yawing  moment  coefficient  with  yaw  rate  due 
to  the  wing  body-combination,  CNRWB,  cannot  be  predicted  from  a 
generalized  method.  However,  the  DATCOM  considers  the  fuselage 
effects  on  this  derivative  as  negligible  and  calculates  the  CNR 
for  wing  alone  to  represent  CNRWB.  In  the  DATCOM  CNRWB  is  given 
as 


CRNWB  =  RCNRCL2  *  CLW2  +  RCNRCDO  *  CDOW  (3.6.2) 

where 


CLW 

RCNRCL2 


RCNRCDO 

CDOW 


is  the  wing  lift  coefficient 

is  the  low  speed  drag-due-to-lift  yaw  damping  parameter 
obtained  from  figure  3.6.1  as  a  function  of  wing  aspect 
ratio  ( ARW ) ,  taper  ratio  (TRW),  sweepback  ( C4SWEEP ) , 
and  mean  aerodynamic  chord  ( MACW ) . 

is  the  low-speed  profile  drag  yaw-damping  parameter 
obtained  from  figure  3.6.2  as  a  function  of  ARW,  C4SWEEP , 
and  MACW 

is  the  wing  profile  drag  coefficient  evaluated  at  the 
appropriate  Mach  number 


3.6.2  CNRV 


The  variation  of  yawing  moment  coefficient  with  yaw  rate  due 
to  the  vertical  tail  is  represented  as 

CNRV  =  -CLALPVO  *  BV  *  SWFV  *  SV/SREF  *  VNO  (3.6.3) 


where 


CLALPVO 

SWFV 

SV 

SREF 

VNO 


is  the  zero  alpha  lift  curve  slope  of  the  vertical  tail 

is  the  sidewash  factor  of  the  vertical  tail 

is  the  area  of  the  vertical  tail 

is  the  reference  area 

is  the  number  of  vertical  tails 


BV  =  2.0  *  (LV  *  C0S( ALPHA ) )2/SPANW2 


(3.6.4) 
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4.0  DYNAMIC  ANALYSIS  ROUTINES 

The  methods  described  below  represent  the  analysis  of  dynamic 
flying  qualities  of  the  aircraft  of  interest.  These  methods  show 
how  to  compute  the  approximate  character istics  for  Short  Period 
Mode,  Phugoid  Mode,  Dutch-Roll  Mode,  and  Rolling  Mode.  The 
dynamic  analysis  methods  presented  here  are  intended  to  give  the 
user  an  approximate,  rigid  body,  non-augmented  dynamic  analysis. 

4.1  SHORT  PERIOD  MODE  EQUATION 

The  short  period  oscillation  of  aircraft  motion  can  be 
characterized  by  the  equation: 

S2  +  2.0  *  ZETA  *  OMEGAN  *  S  +  OMEGAN2  =0  (4.1.1) 

where 

OMEGAN  is  the  undamped  natural  frequency  of  oscillation 
ZETA  is  the  damping  ratio 

TS  is  the  period  of  oscillation 

T12S  is  the  time  to  half  amplitude 

These  characteristics  are  defined  by 

OMEGAN  =  SQRTUMQ  *  LALP/VEL)  -  MALP )  (4.1.2) 

ZETA  =  ((LALP/VEL)  -  MQ  -  MALPD)/(2.0  *  OMEGAN)  (4.1.3) 

TS  =  (2.0  *  PI)/ (OMEGAN  *  SQRT(1.0  -  ZETA2)) 

T12S  =  . 693 1 472/ ( ZETA  *  OMEGAN) 

where 

LALP  =  (G  *  CLALP ) / CLO 

MALP  =  (G  *  CMALP ) / ( MACW  *  CLO  *  (KYS/MACW)2) 

MALPD  =  (G  *  CMALPD )  / ( 2 . 0  *  VEL  *  CLO  *  (KYS/MACW)2) 

MQ  =  (G  *  CMQ ) / ( 2 . 0  *  VEL  *  CLO  *  (KYS/MACW)2) 

and 

G  is  the  force  of  gravity  (ft/sec2) 

VEL  is  the  aircraft  velocity 

CMALD  is  the  variation  of  pitching  moment  coefficient  with  ALPHA 

CMALPD  is  the  variation  of  pitching  moment  coefficient  with  the 

rate  of  change  of  ALPHA 


(4.1.4) 

(4.1.5) 


4 


CLO  is  the  zero  ALPHA  lift  coefficient 

CMQ  is  the  variation  of  pitching  moment  coefficient  with  the 

pitch  rate 

MACW  is  the  mean  aerodynamic  chord  of  the  wing 

KYS  is  the  aircraft  radius  of  gyration  about  the  Y  axis 
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4.2  PHUGOID  MODE  EQUATION 

The  Phugoid  mode  equation  characterizes  a  lower  frequency 
aircraft  oscillation.  The  equation  in  this  case  is  the  same  as 


(v 

A 
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S2  +  2.0  *  ZETA  *  OMEGAN  *  S  +  OMEGAN2  =0  (4.2.1 

Where  the  characteristics  of  the  equation  are  now  represented  by 


OMEGAN 

=  SQRT (2.0)  *  G/VEL 

(4.2.2  ) 

ZETA  = 

SQRT( 2 . 0  )  /  2 . 0  *  CDO/CLO 

(4.2.3) 

TS  =  (2 

.0  *  PI ) / ( OMEGAN  *  SQRT ( 1 . 0  -  ZETA2)) 

(4.2.4  ) 

T 1  2S  = 

.693 1472/ (ZETA  *  OMEGAN) 

(4.2.5  ) 

and 

G 

VEL 

CDO 

CLO 

is  the  force  of  gravity  (ft/sec2) 
is  the  aircraft  velocity  (ft/sec) 
is  the  zero  ALPHA  drag  coefficient 
is  the  2ero  ALPHA  lift  coefficient 
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4.3  DUTCH-ROLL  MODE  EQUATION 

The  Dutch-Roll  mode  equation  is  primary  means  of  analyzing 
aircraft  lateral  stability.  The  characteristic  equation  is  once 
again  given  as; 

S2  +  2.0  *  ZETA  *  OMEGAN  *  S  +  OMEGAN2  =0  (4.3.1) 

with  parameters  defined  as 

OMEGAN  =  SQRT (  ( G  *  CNB ) / ( SPANW  *  CLO  *  ( KZS/SPANW  )  2  )  ) 

ZETA  =  ((G  *  CNR)/ (2.0  *  VEL  *  CLO  *  (KZS/SPANW)2))/ 

(2.0  *  OMEGAN) 

TS  =  (2.0  *  PI)/ (OMEGAN  *  SQRT(1.0  -  ZETA2)) 

T12S  =  . 6 9 3 1 4 7 2 / ( ZETA  *  OMEGAN) 
and 

G  is  the  force  of  gravity  (ft/sec  ) 

VEL  is  the  aircraft  velocity  (ft/sec) 

SPANW  is  the  span  of  the  wing 

KZS  is  the  aircraft  radius  of  gyration  about  the  Z  axis 

CNB  is  the  variation  of  yawing  moment  coefficient  with 

sideslip  angle 

CLO  is  the  aircraft  zero  ALPHA  lift  coefficient 

CNR  is  the  variation  of  yawing  moment  coefficient  with  the 

yaw  rate 


(4.3.2  ) 

(4.3.3  ) 


4.4  ROLLING  MODE  EQUATION 


The  last  dynamic  analysis  equation  characterizes  aircraft 
roll.  This  equation,  known  as  the  Rolling  mode  is  given  by 

P(t)  -  LP  *  BANK(t)  =  LDELA  *  DELTAA 

where 

P(t)  is  the  roll  rate  at  time  (t) 

BANK ( t )  is  the  roll  angle  at  time  (t) 

LP  is  the  roll  damping  parameter 

LDELA  is  the  aileron  power  generating  roll 

DELTAA  is  the  aileron  input 

The  parameters  LP  and  LDELA  are  given  by 

LP  =  (G  *  CLP)/ (2.0  *  VEL  *  CLO  *  ( KXS/SPANW ) 2 ) 

LDELA  =  (G  *  CLDA ) / ( SPANW  *  CLO  *  (KXS/SPANW)2) 

where 

G  is  the  force  of  gravity  (ft/sec^) 

VEL  is  the  aircraft  velocity  (ft/sec) 

KXS  is  the  aircraft  radius  of  gyration  about  the  X  axis 

SPANW  is  the  span  of  the  wing 

CLO  is  the  zero  ALPHA  lift  curve  slope  of  the  aircraft 

CLP  is  the  variation  of  rolling  moment  coefficient  with 

roll  rate 

CLDA  is  the  variation  of  rolling  moment  coefficient  with 

aileron  angle 

From  the  above  results,  the  roll  angle  and  roll  rate  at  a  certain 
time  (t)  after  a  unit  step  input  is 

BANK(t)  =  (-LDELA  *  DELTAA )/LP  *  (t  +  ( 1 /LP )  *  (1.0  -  eLF*t ) ) 

P(t)  =  (-LDELA  *  DELTAA  )/LP  *  (1.0  +  eLP*t ) 

Simultaniously  solving  these  equations  for  (t)  with  a  given 
BANK(t)  and  P(t)  =  0  gives  the  time  required  to  achieve  the 
BANK ( t )  angle.  This  solution  can  only  be  obtained  by  use  of  an 
iterative  method. 


5.0  STATIC  ANALYSIS  METHODS 


The  methods  described  below  allow  for  the  analysis  of 
aircraft  configurations  under  static  conditions.  The 
characteristics  to  be  examined  are  engine  out  conditions,  takeoff 
rotation,  and  crosswind  landing.  The  objective  of  the  methods  is 
to  determine  how  the  aircraft  reacts  in  these  situations. 


5.1  ENGINE  OUT  CONDITIONS 

The  objective  of  this  method  is  to  determine  how  the  aircraft 
performance  degrades  when  one  or  more  of  the  engines  are  out.  The 
affected  systems  are  the  required  rudder  deflection,  the  required 
aileron  deflection,  and  the  minimum  control  velocity. 

First  the  moment  forces  for  the  remaining  engines  (TXM, TXD) 
must  be  determined.  These  are  the  sum  of  the  moments  of  each 
individual  engine.  For  each  engine  the  moments  are  defined  as 

XM  =  THRUST  *  YTAP  (5.1.1) 

XD  =  1.17  *  DQ  *  EIN  *  YTAP  (5.1.2) 

where 

THRUST  is  the  thrust  of  the  engine 

YTAP  is  the  Y  location  of  the  thrust  application  point 

DQ  is  the  atmospheric  dynamic  pressure 

EIN  is  the  engine  inlet  area 

The  Rudder  deflection  required  ( DELRR )  car.  now  be  defined  as 
DELRR  =  ( ( TXM  +  TXD  )  / ( DQ  *  5W  *  SPANW )  ) /CNDELRT  (5.1.3) 

where 

SW  is  the  area  of  the  wing 

SPANW  is  the  span  of  the  wing 

CNDELRT  is  the  variation  of  yawing  moment  coefficient  with 
rudder  angle  (the  sum  of  individual  coefficients  if 
two  tails  are  present) 

The  Aileron  deflection  required  DELAR  is; 

DELAR  =  ( CLBE  *  (TXM  +  TXD ) ) / ( DQ  *  SW  *  SPANW  *  CNB  *  CLDA ) 

(5.1.41 

where 

CLBE  is  the  variation  of  rolling  moment  a  efficient  with 

sideslip  angle 


CNB 


is  the  variation  of  yawing  moment  coefficient  with 
sideslip  angle 

CLDA  is  the  variation  of  rolling  moment  coefficient  with 

aileron  angle 

The  minimum  velocity  for  control  of  the  aircraft  can  be  found  from 
VMIN  =  SQRT ( QMC /  (  .  5  *  DEN))  (5.1.5) 

where 

QMC  =  ABS ( -TXM  -  TXD  )  / ( TN  )  (5.1.6) 

DEN  is  the  atmospheric  density 

TN  =  CLALPVO  *  ALFDR  *  LMACVE  *  DEF  *  SV  (5.1.7) 

ALFDR  is  the  rudder  effectiveness  factor  obtained  from 

f igure  5.1.1 

ALFDR  =  0  for  non-moveable  vertical  tail 

ALFDR  =  1  for  an  all  moveable  vertical  tail 

ALFDR  =  .2923253  +  .4191569  *  CFC  +■  .4816274  *  CFC2  (5.1.8) 

for  a  vertical  tail  rudder 

CFC  is  the  rudder  chord  to  vertical  tail  chord  ratio 

CLALPVO  is  the  vertical  tail  zero  ALFHA  lift  curve  slope 

LMACVE  is  the  mean  aerodynamic  chord  of  the  exposed  vertical 
tail 

DEF  is  the  ta  il  or  rudder  deflection 

SV  is  the  area  of  the  vertical  tail 

If  multiple  vertical  tails  are  present  then  TN  is  summed  for  each 

tail. 


5.2  TAKEOFF  ROTATION 


The  methods  described  here  compute  the  speed  at  which  takeoff 
rotation  occurs  given  the  control  surface  deflections,  and  the 
takeoff  power  setting.  The  first  step  is  the  computation  of  the 
moments  acting  on  the  aircraft.  These  moments  are 

Due  to  aircraft  pitching  moment 

XMAPM  =  CM  *  QTOR  *  SREF  *  MACW  (5.2.1  ) 

where 

CM  is  the  total  aircraft  pitching  moment  coefficient 

SREF  is  the  reference  area 

MACW  is  the  mean  aerodynamic  chord  of  the  wing 

and 

QTOR  =  .5  *  DEN  *  TOVEL2  (5.2.2) 

where 

DEN  is  the  atmospheric  density 

TOVEL  is  the  estimated  takeoff  velocity 

Due  to  engines 

XMEZ  =  -T  *  COS (XI)  *  (-ZTAP  -  ZGLCG ) 

XMEX  =  T  *  SIN(XI)  *  ( XTAP  *  MACW) 

XME  =  XMEX  +  XMEZ  (5.2.3) 

where 

T  is  the  net  thrust  of  the  engine 

XI  is  the  thrust  incidence  angle 

ZTAP  is  the  thrust  application  point  in  the  Z  direction 

XTAP  is  the  thrust  application  point  in  the  X  direction 

ZGLCG  is  the  vertical  distance  between  the  ground  line  and 

the  aircraft  center  line 

These  calculations  are  repeated  for  Body,  Nacelle  type  1 ,  Nacelle 
type  2,  and  lift  engines. 

Due  to  Aircraft  Weight 

XMAW  =  -WEIGHT  *  ( XCG  -  XMG )  *  MACW  (5.2.4) 

where 

XCG  is  the  X  location  of  the  aircraft  center  of  gravity 

XMG  is  the  X  location  of  the  main  landing  gear 
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Due  +-0  wing-body  lift 

LWB  =  (CLOW  *  QTOR  *  SW )  +  ( CLALPWB  *  IW  *  QTOR  *  SREF ) 

XMWBL  =  LWB  *  (XACWB  -  XMG )  *  MACW  (5.2.5) 

where 

LWB  is  the  lift  due  to  the  wing  body 

CLOW  is  the  wing  zero  alpha  lift  curve  slope 

SW  is  the  wing  area 

CLALPWB  is  lift  curve  slope  for  the  wing  and  body  at  ALPHA 

IW  is  the  incidence  of  the  wing 

XACWB  is  the  aerodynamic  center  of  the  wing  body  combination 

Due  to  drag 

CDTAC  =  CDO  +  CDMG  +  CDOTO 
DRAG  =  CDTAC  *  QTOR  *  SREF 

XMAD  =  DRAG  *  ( -ZMACWE  -  ZGLCG )  (5.2.6) 

where 

CDO  is  the  zero  ALPHA  drag  coefficient 

CDMG  is  the  drag  coefficient  due  to  the  main  gear 

CDOTO  is  the  drag  coefficient  for  out  of  ground  effects  and 

high  lift  takeoff  configurations 
ZMACWE  is  the  mean  aerodynamic  chord  of  the  exposed  wing 

Due  to  horizontal  tail  and  elevator 

LHT  =  (CLALPH  *  ( THORR  *  (1.0  -  DEWH  )  -  EWGHT  +  IH  + 

( ADE  *  TELER ) )  *  QTOR  *  NH  *  ASEXH )  +  ( CLOH  *  QTOR  *  ASEXH ) 
XMHT  =  LHT  *  ( XACH  -  XMG)  *  MACW  (5.2.7) 


where 


CLALPH 

THORR 

DEWH 

EWGHT 

IH 

ADE 


is  the 
tail 
>is  the 
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is  the 
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lift  coefficient  due  to  ALPHA  of  the  horizontal 

control  deflection  of  the  horizontal  tail 
downwash  gradient  from  the  wing  to  the  horizontal 

downwash  angle  on  the  horizontal  tail 
incidence  angle  of  the  horizontal  tail 
elevator  effectiveness  parameter  obtained  from 
5.1.1  ( ALFDR ) 


ADE  =  0  for  non-moveable  horizontal  tail 

ADE  =  1  for  an  all  moveable  horizontal  tail 

ADE  =  .2923253  +  .4191569  *  CFC  +  .4816274  *  CFC2  (5.2.8) 

for  a  horizontal  tail  elevator 


TELER  is  the  control  deflection  of  the  elevator 

NH  is  the  dynamic  pressure  ratio  at  the  horizontal  tail 

ASEXH  is  the  exposed  area  of  the  horizontal  tail 


CLOH 


is  the  zero  ALPHA  lift  coefficient  of  the  horizontal 
tail 

XACH  is  the  X  location  of  the  horizontal  tail  aerodynamic 

center 

Due  to  canard  and  canard  flaps 

LC  =  (CLALPC  *  (TCANR  *  (1.0  -  DEWC )  +  IC  +  ( ADCF  *  TCTFR )  )  * 

QTOR  *  NC  *  ASEXC )  +  ( CLOC  *  QTOR  *  ASEXC ) 

XMC  =  LC  *  ( XACC  -  XMG )  *  MACW  (5.2.9) 

where 

CLALPC  is  the  lift  coefficient  due  to  ALPHA  of  the  canard 

TCANR  is  the  control  deflection  of  the  canard 

DEWC  is  the  upwash  gradient  from  the  wing  to  the  canard 

IC  is  the  incidence  of  the  canard 

ADCF  is  the  canard  flaps  effectiveness  factor  obtained  from 

figure  5 . 1 . 1 ( ALFDR  ) 

ADCF  =  0  for  non-moveable  canard 

ADCF  =  1  for  an  all  moveable  canard 

ADCF  =  .2923253  +  .4191569  *  CFC  +  .4816274  *  CFC2  (5.2.10) 

for  a  canard  flap 

TCTFR  is  the  control  deflection  of  the  canard  flap 

NC  is  the  dynamic  pressure  ratio  at  the  canard 

ASEXC  is  the  exposed  area  of  the  canard 

CLOC  is  the  zero  ALPHA  lift  coefficient  of  the  canard 

XACC  is  the  X  location  of  the  canard  aerodynamic  center 

Due  to  acceleration 
R  =  WEIGHT  -  LWB  -  LHT  -  LC 

UDOT  =  ( G  *  ( TTOTAL  -  ( XMEU  *  R)  -  DRAG )) /WEIGHT 

XMA  =  - ( WEIGHT/ G  )  *  UDOT  *  2GLCG  (5-2.11) 

where 

TTOTAL  is  the  total  thrust  of  the  aircraft 

XMEU  is  the  ground  friction  coefficient 

Due  to  High  lift  devices 
LHL  =  ( CLTO  -  CL)  *  QTOR  *  SW 

XMHL  =  -LHL  *  ( XHL  -  (XMG  *  MACW))  (5.2.12) 

where 

CLTO  is  the  lift  coefficient  for  high  lift  take  off 

configurations 

is  the  aircraft  lift  coefficient 
is  the  X  location  of  high  lift  devices 


CL 

XHL 


Now  the  aircraft  total  moment  without  tail(s)  deflection  is 

XMTOTAL  =  XMAPM  +  XMBE  +  XMN1E  +  XMN2E  +  XMLEX  +  XMLEZ  +  XMAW  + 

XMWBL  +  XMAD  +  XMA  +  XMHL  (5.2.13) 


The  above  calculations  are  repeated  for  a  range  of  velocities. 

Two  velocity  points  are  then  located,  the  first  being  the  highest 
velocity  where  the  aircraft  moment  is  greater  than  the  total 
control  moment,  the  second  being  the  lowest  velocity  where  the 
aircraft  moment  is  less  that  the  total  control  moment.  From  these 
two  points  the  takeoff  rotation  is  then  estimated  with 

XT ( 1  )  =  TOVEL ( 1  ) 

XT ( 2 )  =  TOVEL ( 2 ) 

YT1 ( 1 )  =  XMTOTAL ( 1 ) 

YT1 ( 2 )  =  XMTOTAL ( 2 ) 

YT2 ( 1 )  =  XMHT( 1 )  +  XMC( 1 ) 

YT2 ( 2  )  =  XMHT ( 2  )  +  XMHT(2) 

the  velocity  of  takeoff  rotation  ( VTOR )  can  be  computed  by 

TM1  =  ( YT 1 (  1  )  -  YT1 (2)  )  / ( XT ( 1  )  -  XT(2)) 

TM2  =  ( YT2 ( 1 )  -  YT2 ( 2 ) ) / ( XT( 1 )  -  XT(2)) 

B 1  =  YT1  ( 1  )  -  ( TM1  *  XT ( 1  )  ) 

B2  =  YT2 ( 1  )  -  ( TM2  *  XT( 1  )  ) 

and 

TOR  =  ( B2  -  B 1 ) / ( TM 1  -  TM2 )  (5.2.14) 
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5.3  CROSS-WIND  LANDING 


The  methods  described  here  are  used  to  calculate  the  rudder 
deflection  required  to  hold  a  side-slip  angle  of  11.5  degrees  at 
landing  conditions.  The  side-slip  angle  is  found  from  an  assumed 
90  degree  cross-wind  equal  to  .2  times  the  takeoff  velocity  (.2  * 
VTO ) .  The  required  angle  is  found  from 

TAN(B)  =  (VCROSS/V)  =  (.2  *  VTO/VTO)  =  .2  (5.3.1) 

and 

B  =  ARCTAN ( . 2 )  =  11.5  DEGREES 

Now  using  the  equation  of  Equilibrium 

CNB  *  B  =  -CNDELRT  *  DELRCWT  (5.3.2) 

then  the  required  deflection  can  be  found  from 
DELRCWT  =  ( -CNB  *  1  1  . 5  ) /CNDELRT 
where 

CNB  is  the  variation  of  yawing  moment  coefficient  with 

side-slip  angle  evaluated  at  VTO 

CNDELRT  is  the  variation  of  yawing  moment  coefficient  with 
rudder  deflection  evaluated  at  VTO. 


APPENDIX  A  SIGN  CONVENTIONS 


The  sign  convention  used  in  SACP  is  not  the  standard 
aerodynamic  convention.  It  was  changed  to  make  the  conventions 
used  mathematically  correct.  This  change  greatly  simplified  the 
computations  made  in  the  computer  program. 

The  convention  for  moment  arms  is  to  define  negative  arms 
forward  and  above  the  aircraft  with  positive  arms  below  and 
behind.  Figure  A.1  shows  the  arms  in  relation  to  the  aircraft. 
Note  that  there  are  not  arms  depicted  in  the  Y  direction.  The  Y 
axis  is  considered  to  be  positive  in  both  directions  so  components 
on  that  axis  are  not  defined  here. 

The  convention  for  forces  acting  on  the  aircraft  is  that 
positive  forces  act  above  and  behind  the  aircraft,  with  negative 
forces  forward  and  down.  Figure  A. 2  shows  this  relationship. 

The  moment  conventions  are  reversed  from  standard, 
aerodynamic  definitions,  but  in  the  current  configuration  they  are 
mathematically  correct.  With  this  convention,  looking  at  the 
aircraft  from  the  left  wing  tip,  all  positive  moments  act  in  the 
counter  clockwise  direction  and  negative  moments  act  in  the 
clockwise  direction.  Figure  A. 3  details  these  actions  in  each 
quadrant  of  the  axis. 

Figure  A. 4  illustrates  an  example  of  how  the  surface  moment 
arms  are  defined.  The  surface  is  mounted  at  incidence  (is)  with 
forces  of  lift(L)  and  drag  (D)  acting  upon  it.  In  this  case  the 
moment  arms  are  defined  to  be 

X  =  (XACS  *  COS  B )  +  (ZACS  *  SIN  g) 

Z  =  ( XACS  *  SIN  B )  -  ( ZAC  s  *  COS  g) 
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APPENDIX  B  GEOMETRY 


Though  not  an  aerodynamic  method,  the  proper  calculation  of 
geometry  is  essential  to  the  accurate  development  of  aerodynamic 
equations.  The  following  section  describes  the  input  geometry  for 
the  wing,  and  the  subsequent  geometry  calculations .  The  results 
obtained  here  can  be  extended  directly  to  the  other  aircraft 
surfaces . 

Figures  1  show  the  wing  geometry  inputs  which  are  defined  as 

follows : 

NATW  the  wing  airfoil  option  flag 

NPANSW  the  number  of  wing  panels  (4  max) 

CSAW  the  average  body  cross-section  area  along  the  wing 

root  chord  (SQFT) 

BDW  the  body  diameter  along  the  wing  root  chord  (ft) 

XTCMW  the  location  of  maximum  t/c  along  the  wing 

mean  aerodynamic  chord  (%  MAC) 

IW  the  wing  incidence  (DEG) 

IFLEXW  the  wing  leading  edge  extension  bluntness  flag 

SEXTW  the  wing  leading  edge  extension  planform 

area  (SQFT) 

SLEXTW  the  wing  leading  edge  extension  leading 

edge  sweep  ( DEG  ) 

XEXTW  the  wing  leading  edge  extension  centroid  x 

station  ( ft ) 

CREW  the  wing  exposed  root  chord  length  ( ft ) 

TCREW  the  wing  T/C  of  exposed  root  chord  leading 

edge  ( % ) 

XCREW  the  wing  x  location  of  exposed  root  chord 

leading  edge  ( ft ) 

YCREW  the  wing  y  location  of  exposed  root  chord 

leading  edge  (ft) 

ZCREW  the  wing  z  location  of  exposed  root  chord 

leading  edge  (ft) 

WSPLEP(NP)  the  leading  edge  sweep  of  each  panel  of  the  wing 
( DEG ) 

WTOCP(NP)  the  t/c  at  the  tip  chord  of  each  panel  of  the  wing 
( %  ) 

WSPANP(NP)  the  span  of  each  panel  of  the  wing  (ft) 

WTCP(NP)  the  tip  chord  length  of  each  panel  of  the 
wing  ( ft ) 

WOFFP(NP)  the  leading  edge  offset  length  (-  forward)  of  the 
wing  ( ft ) 

WCAMP(NP)  the  design  lift  camber  coefficient  of  each 
panel  of  the  wing 

WTWISTP(NP)  the  twist  of  each  panel  (-  for  washout)  of  the 
wing 

WDIHP(NP)  the  dihedral  angle  of  each  panel  (+  up)  of  the 
wing  (DEG) 

FIGURE  WING  PANEL  INPUTS 


With  these  inputs  some  additional  properties  can  be  computed  by 


SXXW(NP)  the  exposed  area  of  each  panel 
SXXW(NP)  =  (WCTP(NP)  +  WCTP(N-I))  *  WSPANP(NP) 

WCTP ( 0  )  =  CREW 

ASEXW  the  exposed  area  of  the  wing 

ASESW  =  SXXN(NP) 

ARWPE(NP)  the  exposed  aspect  ratio  of  each  panel 
ARWPE(NP)  =  (2.  *  WSPANP ( NP )  ) 2 /SXXW ( NP  ) 

TRWPE(NP)  the  exposed  aspect  ratio  of  each  panel 
TRWPE(NP)  =  WCTP (NP) /WCTP ( NP- 1 ) 

XCRW(NP)  the  exposed  root  chord  of  each  panel 
XCRW(NP)  =  XCRS ( NP- 1 )  +  TAN( WSPLEP ( NP- 1 ) )  *  WSPANP(NP-l)  + 
SOFFP(NP) 

XCRW('i)  =  XCREW 

SPANWE  the  exposed  wing  span 
SPANWE  =  WSPANP (NP)  *  2. 

SPANW  the  wing  span 
SPANW  =  SPANWE  +  YCREW  *  2. 

WSPC4P( NP )  the  quarter  chord  sweep  of  each  panel 
WSPC4P(NP)  =  ATAN( ( ( ( WSPANP (NP)  *  2.  *  (TAN  ( WSPLEP ( NP ))  )  + 

(.5  *  WCTP ( NP ) )  -  (.5  *  WCTP (NP-1 ))/( WSPANP  *  2.) 

WSPC2P(NP)  the  half  chord  sweep  of  each  panel 

WSPC2P( NP )  =  ATAN(  (( (WSPANP ( NP ) *  2.  *  (TAN  ( WSPLEP ( NP  ))  )  + 

( WCTP ( NP )  )  -  ( WCTP ( NP  -!))/( WSPANP  *  2.) 

WSPTEP ( NP )  the  trailing  edge  sweep  of  each  panel 
WSPTEP(NP)  =  ATAN(  (  (  ( WSPANP (NP)  *  2.  *  ( TAN ( WSPLEP ( N2  ))  )  + 

(2.  *  WCTP(NP))  -  (2.  *  WCTP ( NP  -!))/( WSPANP  *  2.) 

WSPMTP(NP)  the  maximum  thickness  sweep  of  each  panel 
WSPMTP ( NP )  =  AT AN( (( (WSPANP (NP)  *  2.  *  ( TAN ! WSPLEP ( N2  ))  )  + 

( XTCMW  *  2.  *  WCTP ( NP )  )  -  (XTCMW  *  2  *  WCTP(N?  -  1)) 
/(WSPANP  *  2) 

SUM  1  =  TAN ( WSPLEP ( NP )  )  *  SXXW(NP) 

SUM2  =  TAN( WSPC4P(NP ) )  *  SXXW(NP) 

SUM3  =  TAN( WSPC2P( NP ) )  *  SXXW(NP) 

SUM4  =  TAN ( WSPTEP ( NP ) )  *  SXXW(NP) 

SUM5  =  TAN( WSPMTP ( NP ) )  *  SXXW(NP) 

SPLEWR  the  total  wing  leading  edge  sweep 
SPLEWR  =  ATAN( SUM1 /ASEXW) 
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SPC4WR  the  total  wing  quarter  chord  sweep 
SPCRWR  =  ATAN ( SUM2 / ASEXW ) 

SPC2WR  the  total  wing  half  chord  sweep 
SPC2WR  =  ATAN ( SUM 3 / ASEXW ) 

SPTEWR  the  total  wing  trailing  edge  sweep 
SPTEWR  =  ATAN ( SUM4 / ASEX ) 

SPMTWR  the  total  wing  maximum  thickness  sweep 
SPMTWR  =  ATAN ( SUM5/ ASEX ) 

WMACP(NP)  the  mean  aerodynamic  chord  of  each  panel 
WMACP(NP)  -2/3  *  ( WCTP ( NP  -  1)  -  WOFFP(NP))  *  (1.  +  TRWPE(NP) 

+  TRWPE( NP ) 2 ) / (  1  .  +  TRWPE ( NP  )  ) 

AWET(NP)  the  wetted  area  of  each  panel 

AWET(NP)  =  SXXW(NP)  *  (2  +  .1843  *  WTOCP(NP)  +  1.5268  * 

STOCP ( NP  ) 2  -  .8395  *  STOCP(NP)3) 

ARW  the  aspect  ratio  of  the  wing 
ARW  =  SPANW2/SW 

ARWE  the  exposed  aspect  ratio  of  the  wing 
ARWE  =  SPANWE2 /ASEXW 

XCR23W  the  X  location  of  the  2/3  chord  point 
XCR23W  =  XCREW  +  (2/3  *  CREW) 

ZCGCRW  the  Z  distance  between  the  center  of  gravity  and  the  root 
chord  of  the  wing 

ZCGCRW  =  ZCREW  -  ((.25  *  CREW)  *  SIN(IW)) 
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